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The  ability  to  accurately  predict  temperatures  in  a  circulating  wellbore  is 
blv  more  critical  today  than  at  any  time  in  the  history  of  drilling  oil  wells, 
o  the  increasing  depths  to  which  wells  are  being  drilled,  accurate  methods  and 
dures  are  needed  to  predict  fluid  temperatures  during  circulation  operations, 
methods  and  procedures  would  be  used  to  provide  accurate,  data  for  designing 
ing  mud  systems  and  cement  jobs  for  these  d  "p ,  hot.  wells. 

Though  the  problem  of  accurately  predicting  temperatures  in  a  wellbore  has 
ed  since  engineers  first  began  searching  for  oil  at  greater  depths,  relatively 
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Numerous  large-scale  temperature  simulators  have  been  developed  by  many 
researchers.  However,  these  simulators  are  restricted  in  their  applicability  because 
of  their  complexity  and  limited  portability.  As  an  alternative,  a  model  has  been 
developed  for  the  personal  computer  which  is  portable,  flexible,  and  easy  to  use 

The  model  logically  organizes  the  required  data  into  categories.  Previous  tfafei 
files  can  be  re-entered  and  modified  as  need  be.  Menus  are  used  to  control  the 
program  and  select  different  options.  A  liberal  portion  of  on-screen  computer 
graphics  is  included  to  insure  that  the  model  is  truly  "user  friendly."  The  results 
of  the  simulator  have  been  compared  with  published  information  from  other  models. 

The  comparison  shows  that  a  small  scale  model  can  be  designed  for  and  operated  on  a 
personal  computer,  include  a  considerable  amount  of  user  interactiveness,  and 
maintain  realism  and  accuracy. 

Also,  the  results  of  a  parametric  sensitivity  analysis  carried  out  using  the 
model  are  discussed.  The  outcome  demonstrates  that  certain  parameters  which  may  have 
been  ignored  in  previous  work,  have  a  significant  effect  on  wellbore  temperature 
profiles . 
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ABSTRACT 


A  PC  Simulation  of  Heat  Transfer  and  Temperature 
Distribution  in  a  Circulating  Wellbore 
(December  1987) 

Robert  Duane  Pierce,  B.S.,  South  Dakota  School 
of  Mines  and  Technology 

Chair  of  Advisory  Committee:  Dr.  Hans  C.  Juvkam-Wold 

The  ability  to  accurately  predict  temperatures  in  a  circulating 
wellbore  is  probably  more  critical  today  than  at  any  time  in  the 
history  of  drilling  oil  wells.  Due  to  the  increasing  depths  to  which 
wells  are  being  drilled,  accurate  methods  and  procedures  are  needed  to 
predict  fluid  temperatures  during  circulation  operations^.  These 
methods  and  procedures  would  be  used  to  provide  accurate  data  for 
designing  drilling  mud  systems  and  cement  jobs  for  these  deep,  hot 
wells . 

Though  the  problem  of  accurately  predicting  temperatures  in  a 
wellbore  has  existed  since  engineers  first  began  searching  for  oil  at 
greater  depths,  relatively  little  has  been  done  to  address  this 
problem.  In  the  past  it  has  been  convenient  to  ignore  temperature 
gradients,  largely  because  no  practical  means  for  estimating  wellbore 
temperature  profiles  have  been  available'1.  Determining  accurate 
temperature  distributions  in  a  circulating  wellbore  is  very  important 
for  many  aspects  of  drilling,  completion,  production,  and  injection. 

Numerous  large-scale  temperature  simulators  have  been  developed 


by  many  researchers.  However,  these  simulators  are  restricted  in 
their  applicability  because  of  their  complexity  and  limited 
portability.  As  an  alternative,  a  model  has  been  developed  for  the 
personal  computer  which  is  portable,  flexible,  and  easy  to  use. 

The  model  logically  organizes  the  required  data  into  categories. 
Previous  data  files  can  be  re-entered  and  modified  as  need  be.  Menus 
are  used  to  control  the  program  and  select  different  options.  A 
liberal  portion  of  on-screen  computer  graphics  is  included  to  insure 
that  the  model  is  truly  "user  friendly."  The  results  of  the  simulator 
have  been  compared  with  published  information  from  other  models.  The 
comparison  shows  that  a  small  scale  model  can  be  designed  for  and 
operated  on  a  personal  computer,  include  a  considerable  amount  of  user 
interactiveness,  and  maintain  realism  and  accuracy. 

Also,  the  results  of  a  parametric  sensitivity  analysis  carried 
out  using  the  model  are  discussed.  The  outcome  demonstrates  that 
certain  parameters  which  may  have  been  ignored  in  previous  work,  have 
a  significant  effect  on  wellbore  temperature  profiles. 
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The  objective  of  this  research  is  to  develop  a  small  scale 


computer  program  for  use  on  a  standard  personal  computer  which  will 


accurately  model  conductive  and  convective  heat  transfer  processes  in 


a  flowing  wellbore.  Ultimately,  this  model  will  assist  in  providing  a 


greater  degree  of  accuracy  in  the  prediction  of  temperature  in  a 


flowing  wellbore .  The  computer  model  will  incorporate  only  those 


parameters  that  have  a  significant  effect  on  the  system  dynamics, 


maintaining  a  high  degree  of  accuracy  and  realism.  The  proposed 


computer  model  will  be  very  interactive  and  informative,  using  on¬ 


screen  menus  and  graphics .  A  flow  chart  of  the  model  is  found  at 


Fig.l.  The  resulting  simulator  will  be  relatively  fast  running  and 


flexible,  making  it  an  effective  tool  for  well  planning  or  for 


designing  various  cement  jobs.  The  program  will  be  written  on  a  5.25 


inch  floppy  disk,  making  it  usable  in  the  office  or  at  the  well  site. 


The  ability  to  accurately  predict  temperatures  in  a  circulating 


wellbore  is  probably  more  critical  today  than  at  any  time  in  the 


history  of  drilling  oil  wells.  Due  to  the  increasing  depths  to  which 


wells  are  being  drilled,  accurate  methods  and  procedures  are  needed  to 


predict  fluid  temperatures  during  circulation  operations. ^  These 


methods  and  procedures  would  be  used  to  provide  adequate  data  for 


correctly  designing  drilling  mud  systems  and  cement  jobs  for  these 


deep,  hot  wells.  Though  the  problem  of  accurately  predicting 


This  thesis  follows  the  style  of  the  Journal  of  Petroleum  Technolot 
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Fig.  1  ■  Computer  program  flow  chart. 
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temperatures  in  a  wellbore  has  existed  since  engineers  first  began 
searching  for  oil  at  greater  depths,  relatively  little  has  been  done 
to  address  this  problem.  The  present  trend  of  drilling  deeper  and 
more  expensive  wells  requires  a  more  accurate  knowledge  of  the 
variables  involved  in  the  operation,  such  as  temperature.  In  the  past 
it  has  been  convenient  to  ignore  temperature  gradients,  largely 
because  no  practical  means  for  estimating  wellbore  temperature 
profiles  have  been  available. ^  Determining  accurate  temperature 
distributions  in  a  circulating  wellbore  is  very  important  for  many 
aspects  of  drilling,  completion,  production,  and  injection.  A  few  of 
the  applications  that  require  an  understanding  of  the  down-hole 
temperature  history  in  a  well  include:3 

1.  Cement  composition,  placement,  and  setting  time. 

2.  Drilling  mud  and  annulus  fluid  formulation. 

3.  Packer  design  and  selection. 

4.  Logging  tool  design  and  log  interpretation. 

5.  Wax  deposition  in  production  tubing. 

6.  Corrosion  in  tubing  and  casing. 

7.  Thermal  stresses  in  casing  and  tubing. 

8.  Permafrost  thawing  and  refreezing. 

9.  Wellhead  and  production  equipment  design. 

10.  Drill  bit  design. 

11.  Elastomer  and  seal  selection. 

Of  course,  many  other  possible  applications  may  exist.  Two 
interesting  possibilities  for  computer  modeling  are  computing  the 


undisturbed  formation  temperatures  from  flowing  temperature  stream 
measurements  and  predicting  abnormal  pressure  zones  from  fluid 
temperature  changes  while  drilling.  Analytical  means  must  be  used  to 
model  the  system  in  question  since  temperature  data  can  only  be 
obtained  by  physical  measurement  at  isolated  points  within  the 
wellbore.  Such  methods  have  formed  the  basis  for  previous  work  in 
this  area. 

As  early  as  1941,  scientists  and  engineers  have  studied  the 
effect  of  temperature  on  drilling  fluids  and  cements. It  has  been 
well  documented  in  the  literature  that  downhole  temperatures  in  a 
wellbore  which  are  recorded  during  routine  logging  operations,  do  not 
measure  true,  static  formation  temperature.  Due  to  the  cooling 
effects  of  fluid  circulation,  recorded  temperatures  can  be  25°F.  to 
70°F.  lower  than  true,  static  formation  temperature . ^  Many  models  and 
methods,  both  theoretical  and  empirical,  designed  to  approximate 
temperature  profiles  in  the  wellbore  under  various  conditions,  have 
been  proposed  since  the  late  1950' s.^'^l  These  methods  range  from 
very  simple  mathematical  relationships  and  empirical  charts,  which  can 
be  performed  by  hand  calculations,  to  large,  complex,  finite- 
difference  simulators,  suitable  for  use  on  high-speed,  main-frame 
digital  computers.  However,  very  little  recent  data  has  been  gathered 
to  verify  any  of  these  theoretical  models . ^  Also,  some  work  has  been 
done  to  evaluate  the  effects  of  temperature  on  subsurface  equipment  in 
the  wellbore . ^2 - 28  Despite  the  advances  in  reservoir  simulation,  the 
real  use  and  application  of  presently  available  models  has  been 


limited  to  facilities  with  large,  main- frame  computers,  such  as 
training  centers. 

In  1967,  Annis^  published  the  results  of  an  extensive  study  of 
the  effects  of  high  temperature  on  the  flow  properties  of  water-base 
drilling  fluids.  His  investigation  of  the  influence  of  time, 
temperatures  up  to  300°F,  and  mud  composition  on  critical  drilling 
fluid  properties  found  that  gel  strengths  tended  to  be  more  sensitive 
than  viscosities  to  changes  in  temperature  and  mud  composition.  The 
gel  strength  of  bentonitic  clay  drilling  fluids  at  high  temperature 
was  not  predictable  from  measurements  made  at  low  temperatures. 

When  an  engineer  is  designing  a  drilling  mud  program  or  a  cement 
slurry  for  setting  casing,  many  factors  must  be  considered,  such  as 
the  density  of  the  fluid,  fluid  loss  control,  fluid  viscosity,  and  the 
deterioration  of  the  fluid  from  encountering  high  temperatures.  In 
the  case  of  designing  a  cement  job,  pumping  time  and  compressive 
strength  must  also  be  considered.  Depending  on  the  conditions 
encountered  in  individual  wells,  other  factors  may  have  to  be 
considered. 

Pumping  time  is  a  primary  factor  in  obtaining  a  good  cement  job 
and,  as  wells  are  drilled  to  greater  depths,  this  property  becomes 
even  more  important.  This  is  because  we  are  experiencing  ever 
increasing  bottomhole  temperatures.  To  design  a  good  cement  job,  we 
must  have  the  most  accurate  temperature  data  possible.  We  must  design 
for  sufficient  pumping  time  for  the  cement  slurry  to  be  safely  placed 
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in  the  well,  but  the  slurry  cannot  be  overly  retarded  as  this  will 
hinder  the  development  of  satisfactory  compressive  strength. 

The  prediction  of  wellbore  temperature  during  circulation  cannot 
be  accomplished  easily.  This  is  because  it  is  a  complex  function  of 
wellbore  geometry,  geothermal  gradient,  circulation  rate  and  time, 
fluid  properties,  and  film  heat  transfer  coefficients.^  These 
temperature  estimates  which  we  are  seeking  to  determine  during 
drilling  and  circulation  require  an  understanding  of  the  thermal 
processes  which  occur  in  the  wellbore.  The  information  gathered  may 
be  used  to  predict  the  behavior  of  cement  slurries,  drilling  muds,  and 
drilling  equipment  in  the  severe  temperatures  of  the  downhole 
environment . 

A  further  aim  of  this  research  was  to  use  the  computer  simulator 
to  simulate  wellbore  temperature  profiles  so  that  a  thorough 
parametric  sensitivity  analysis  could  be  carried  out  to  determine 
which  parameters  have  the  most  significant  effect  on  the  wellbore 
temperature . 

STATEMENT  OF  THE  PROBLEM 

Although  the  simulation  of  temperature  distributions  in  a 
circulating  wellbore  has  many  fundamental  uses  in  modern  drilling 
fluid  and  equipment  design,  simulation  is  poorly  understood  by 
engineers  in  the  field.  This  is  attributable  to  the  present  nature  of 
heat  and  temperature  simulation  methods  and  procedures.  The  major 
disadvantage  to  most  temperature  prediction  routines  is  that,  for 
example,  some  of  the  methods  currently  available  in  the  literature  are 


based  on  data  and  assumptions  which  do  not  approximate  individual 
wells.  Most  of  the  digital  computer  simulat.  i  s  requite  a  main  frame 
computer  to  operate,  and  are  not  portable  to  the  drill  site 
Therefore,  an  engineer  on  site  may  have  temperature  predictions 
available  to  him  based  on  well  design  factors,  but  has  no  convenient 
means  of  obtaining  real-time  data  on  actual  drilling,  production,  or 
injection  operations. 

Although  these  large-scale  models  are  generally  accurate  and 
effective,  a  more  flexible  approach  to  temperature  distribution 
simulation  is  needed.  For  a  simulator  to  have  broad  applications,  it 
should  be  portable,  accurate,  and  easy  to  use.  Highly  portable 
simulators  would  permit  a  model  to  be  used  anywhere  from  the  classroom 
to  the  office  to  the  drilling  rig  floor.  Accuracy  and  realism  is 
essential  to  allow  simulation  of  various  drilling  and  producing 
conditions.  This  should  be  accomplished  while  producing 
representative  results.  Lastly,  a  simulator  is  not  a  genuinely  useful 
device  unless  it  is  easy  to  use  by  those  who  have  a  need  for  its 
application.  The  requirements  of  the  model  should  be  clear  and  well 
defined,  and  the  results  from  the  simulation  should  be  easy  to 
understand  and  evaluate. 

REVIEW  OF  CURRENT  LITERATURE 

In  1941,  Farris^*  presented  what  is  probably  the  earliest  study  of 
wellbore  temperature.  In  his  paper,  he  developed  charts  and  special 
procedures  for  predicting  the  bottomhole  temperatures  of  circulating 
wells  based  on  measurements  in  five  shallow  Gulf  Coast  wells.  Major 


weaknesses  associated  with  these  charts  and  procedures  prompted 
further  research  into  a  more  precise  method  for  estimating  circulating 
temperatures.  It  was  not  until  some  thirteen  years  later  that  more 
detailed  studies  of  cement  strength  and  setting  time  were  conducted. 

In  March  1954,  Swayze®  presented  the  results  of  a  committee  study  on 
the  effects  of  high  pressures  and  temperatures  on  the  strength  of  oil 
well  cements.  In  April,  O'Neal  and  Benischek®  presented  their  study 
of  the  setting  time  of  cements  as  it  is  affected  by  high  temperatures 
and  pressures. 

In  a  classic  paper,  Ramey®  presented  an  approximate  solution  to 
the  wellbore  heat  transmission  problem  involved  in  the  injection  of 
hot  or  cold  fluids  into  injection  wells.  This  analysis  allowed  for 
the  calculation  of  the  temperature  distribution  in  these  injection 
wells.  His  solution  permitted  estimates  of  the  temperature  of  fluids, 
tubing  and  casing  as  a  function  of  time  and  depth.  Later,  he 
presented  a  paper®®  which  expanded  on  this  derivation  to  give  the  rate 
of  heat  loss  from  the  well  to  the  formation. 

Edwardson  et  al .  ^  developed  a  method  which  has  its  basis  in  the 
mathematical  solution  of  the  differential  equation  of  heat  conduction. 
The  solution  of  this  equation  is  presented  in  a  series  of  graphs  which 
are  used  to  determine  formation  temperature  disturbances  at  various 
radii  for  arbitrary  mud  circulation  conditions.  This  work  was  more 
concerned  with  temperature  distribution  in  the  formation  as  a  result 
of  drilling  fluid  circulation,  and  does  not  allow  for  the  direct 


calculation  of  the  wellbore  temperature  profile. 
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Tragesser  et  al .  u  presented  a  calculation  technique  which 
provided  circulating  temperatures,  as  functions  of  time,  at  varying 
depths  in  both  the  casing  and  annulus.  This  was  an  expansion  of  the 
method  of  Edwardson  et  al .  ^  so  as  to  allow  for  such  variables  as 
depth,  pumping  rate,  hole  diameter,  and  so  forth.  However,  both  of 
these  techniques  are  based  on  an  assumed  formation  temperature 
distribution  after  circulation  and  provided  no  means  of  calculating 
the  wellbore  temperature  distribution  directly. 

Holmes  and  Swift^  obtained  a  steady-state  solution  to  the 
wellbore  heat  transfer  equation  for  the  heat  transfer  between  the 
fluids  in  the  annulus  and  the  fluids  in  the  drill  pipe.  This  was 
combined  with  an  approximate  equation  for  the  transient  heat  transfer 
between  the  fluid  in  the  annulus  and  the  formation.  The  authors  state 
that  the  approximate  method  is  adequate  since  the  total  heat  transfer 
between  the  two  fluids  is  much  greater  than  between  the  annulus  fluid 
and  the  formation.  The  low  heat  transfer  between  the  annulus  fluid 
and  the  formation  is  a  result  of  the  relatively  low  thermal 
conductivity  of  the  formation  and  the  film  resistance  to  heat  transfer 
formed  at  the  interface  of  the  mud  and  the  rock.  Temperatures  can  be 
calculated  as  a  function  of  well  depth,  mud  circulation  rate, 
circulating  fluid  attributes,  reservoir  characteristics,  and  wellbore 
and  drill  pipe  size.  It  is  the  work  of  Holmes  and  Swift^,  which  is 
simple  and  readily  adaptable  to  this  study,  that  is  the  basis  of  the 
research  described  in  this  thesis. 
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Raymond^  developed  a  method  for  predicting  temperature 
distributions  for  both  transient  and  pseudosteady- state  conditions. 

He  developed  a  generalized  technique  for  calculating  the  entire 
temperature  distribution  in  the  system  and  a  general  method  for 
predicting  bottomhole  fluid  temperatures  during  circulation.  Raymond 
advanced  the  use  of  the  principle  of  superposition  and  the  Hurst  and 
van  Everdingen  functions  to  solve  numerically  for  unsteady- state 
conditions.  He  insisted  that  the  pseudosteady-state  solution  was 
adequate  for  all  applications . ^ 

In  1977,  Fertl  and  Wichmann^  introduced  a  simple  and  rapidly 
applied  technique  for  analyzing  maximum  bottomhole  temperatures,  which 
are  recorded  during  well  logging  operations,  to  determine  static 
formation  temperature.  The  method  requires  the  use  of  a  maximum 
recording  thermometer  on  each  logging  run,  plus  information  concerning 
circulating  time  and  time  that  the  logging  tool  was  last  on  the  bottom 
of  the  borehole. 

Most  recent  research  has  dealt  with  means  of  solving  unsteady 
state  equations  formulated  by  Raymond,  or  variations  on  these,  using 
finite  difference  techniques . ^ ^ • 15-21  Although  the  earlier  methods 
are  simple  and  easy  to  apply,  they  are  not  particularly  accurate.  The 
more  recent  computer  models  involve  solving  the  finite  difference 
equations  describing  the  heat  transfer  by  iterative  methods.  Such 
methods  have  the  obvious  disadvantages  of  long  solution  times  and  the 
simultaneous  problems  of  stability  and  accuracy. 


In  1968,  Keller  and  Couch^  studied  the  damage  to  production 
wells  caused  by  the  high  temperatures  associated  with  in- situ 
combustion  projects.  Injection  of  water  down  the  annulus  of  hot  wells 
has  been  successfully  utilized  for  cooling  and  prolonging  the  life  of 
such  wells.  They  developed  a  mathematical  model  that  was  acceptable 
for  predicting  the  effectiveness  of  cooling  water  in  a  particular 
application. 

Keller  ec  al.^  developed  a  model  which  described  the  two- 
dimensional  transient  heat  transfer  in  and  around  a  wellbore.  Their 
results  demonstrated  that  the  use  of  steady- state  solutions  previously 
published  gave  good  estimates  of  the  circulating  mud  temperatures. 

The  transient  solution  that  was  presented  in  their  paper  was  more 
suited  for  matching  temperature  logs. 

Sump  and  Williams^  presented  an  improved  model  which  they  stated 
was  more  accurate  than  the  use  of  API  correlations  or  other  numerical 
techniques.  Their  model  was  based  on  an  existing  mathematical  model 
of  the  wellbore  and  the  formation  near  the  wellbore.  This  existing 
model  was  improved  by  altering  the  film  heat  transfer  coefficients  and 
the  formation  thermal  conductivity.  The  coefficients  were  modified  to 
simultaneously  minimize  the  difference  between  measured  and  predicted 
temperatures  for  six  wells. 

As  computer  models  have  become  more  sophisticated,  numerous 
authors  have  presented  simulators  which  model  even  more  complex  well 
conditions.  An  excellent  example  of  this  is  the  study  by  Oster  and 
Scheffler^  which  described  their  model  for  determining  the 


temperature  distribution  in  a  circulating  drilling  fluid  when  aquifers 
are  present  in  the  formation.  The  depth  of  the  aquifer  relative  to 
the  well  depth  was  shown  to  be  an  important  parameter. 

As  computers  have  become  faster  and  more  able  to  solve 
increasingly  complex  problems,  temperature  distribution  simulators 
have  become  more  complicated.  In  1980,  Wooley-^  presented  a  complex 
model  for  predicting  downhole  wellbore  temperatures  in  flowing  of 
shut-in  fluid  streams,  in  casing  and  cement,  and  in  formations. 

Flowing  options  include  injection/production,  forward/reverse 
circulation,  and  drilling. 

In  1981,  Cline^  formulated  a  mathematical  model  of  the 
temperature  distribution  inside  a  circulating  wellbore.  He  assumed 
that  the  heat  transfer  processes  between  the  formation  and  the 
circulating  fluid  in  the  wellbore  occur  in  a  pseudo-steady  manner.  A 
closed  form  solution  to  the  counterflow  energy  equations  was  developed 
in  a  manner  which  allows  for  arbitrary  temperature  distributions  in 
the  formation. 

Marshall  and  Bentsen^  developed  a  computer  model  which  uses  a 
direct  solution  technique  to  solve  the  finite  difference  equations 
describing  transient  heat  transfer  in  the  wellbore.  The  authors 
claim  that  their  solution  technique  is  considerably  more  efficient 
than  those  used  in  earlier  studies. 

In  1984,  Duda^  presented  computer  calculations  of  wellbore 
transient  temperatures,  using  the  geothermal  wellbore  thermal 
simulator  code  GE0TEMP2 .  These  calculations  were  made  on  four  well 
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model's.  Also  In  1984,  Corre  ec  al  developed  a  computer  model  which 
could  be  used  for  determining  temperature  profiles  in  and  around  a 
wellbore  during  drilling  phases.  This  model  has  been  used  to  design 
mud  and  drilling  programs  to  ensure  better  hole  stability. 

While  various  authors  were  working  to  more  accurately  model  and 
predict  downhole  circulating  mud  temperature  distributions,  other 
authors  were  studying  the  effects  of  temperature  on  subsurface 
equipment . 

Leutwyler  and  Bigelow^  made  a  detailed  study  of  the  effects  of 
temperature  change  on  tubing,  casing,  and  associated  downhole 
equipment  as  a  result  of  injection  or  production  of  fluids  at 
temperatures  substantially  higher  than  the  surrounding  formations. 

They  presented  a  method  for  making  an  approximation  of  the  casing 
temperature  and  for  analyzing  the  buckling  criteria. 

In  1966,  Leutwyler  provided  an  excellent  review  of  much  of  the 
information  leading  to  the  prediction  of  heat  losses  from  tubing 
string  to  the  wellbore.  He  also  reviewed  the  line  source  solution  to 
the  diffusivity  equation  and  examined  some  of  the  application  criteria 
in  light  of  the  unsteady - s tate  conditions. 

Dase  and  Heyt^  studied  the  thermal  protection  of  wellbore 
casings.  Their  analytical  study  of  heat  transfer  through  the  wellbore 
during  hot  fluid  stimulation  was  to  determine  the  effect  of  a 
concentric  radiation  shield  and  forced  convection  gas  flow  as  methods 
for  the  reduction  of  casing  temperature  and  associated  thermal  stress 
which  may  lead  to  failure. 
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DISCUSSION  OF  THE  MODEL 


As  stated  earlier,  the  primary  objective  of  this  research  was  to 
design  a  totally  "user  friendly"  computer  program  to  assist  the 
drilling  engineer  or  the  production  engineer  to  more  accurately  model 
the  temperature  and  heat  transfer  processes  that  are  occurring  in  the 
circulating  wellbore.  During  the  course  of  developing  this  program, 
the  author  was  able  to  include  pressure  loss  calculations.  The 
simulator  should  be  flexible,  easy  to  use,  provide  simple,  effective 
output  which  is  easy  to  use ,  and  require  a  minimum  amount  of  data  and 
time  to  operate.  However,  the  quality  and  accuracy  of  the  results 
should  not  be  compromised.  To  make  the  model  flexible,  it  is  designed 
for  use  on  a  personal  computer.  This  operating  environment  encourages 
the  use  of  floppy  disks  and  minimal  program  size.  The  key  advantage 
to  the  use  of  floppy  disks  is  their  portability.  Because  of  the 
capacity  of  floppy  disks  and  the  variable  memory  sizes  of  personal 
computers,  the  simulator  in  question  had  to  have  reasonable  memory 
requirements.  The  final  version  of  this  heat  transfer  and  temperature 
gradient  simulator  has  a  memory  requirement  of  less  than  50  kilobytes 
and  can  be  operated  on  a  machine  having  256  kilobytes  of  RAM  memory. 

To  enhance  the  practical  nature  of  the  model,  it  was  vital  to 
create  a  set  of  subroutines  that  acted  as  a  complete  package.  This 
suggests  that  the  simulator  be  "user  friendly",  i.e.,  interactive, 
self  explanatory,  flexible,  and  easy  to  use.  To  realize  these 
objectives,  it  was  necessary  to  employ  a  liberal  amount  of  on-screen 
graphics  and  menus.  The  program  HEAT. BAS  creates  the  data  file  which 


describes  the  circulating  well  and  the  heat  and  temperature  processes 
within  it.  The  data  are  logically  separated  by  category  and  are 
presented  as  "pages"  of  information.  Each  "page"  describes  the 
variables,  their  appropriate  units,  and  the  default  values  that  are 
built  into  the  main  program.  Also,  each  "page"  has  the  current  date 
and  time  available  for  the  user.  Menus  let  the  program  user  move 
quickly  between  these  "pages". 

In  order  to  use  this  program,  the  personal  computer  should  also 
have  the  compiler  QUICKBASIC.EXE  in  the  RAM.  Once  this  is 
accomplished,  the  program  HEAT. BAS  can  be  loaded  into  QUICK. BASIC  and 
the  heat  transfer  program  can  be  run.  Upon  executing  HEAT. BAS,  the 
first  "page"  of  information  that  appears  on  the  terminal  screen  is  the 
title  block  with  the  name  of  the  program  and  the  name  of  the  author, 
as  illustrated  in  Fig.  2.  The  user  is  prompted  to  press  any  key  to 
move  to  the  next  "page" . 

Having  pressed  a  key,  the  second  "page"  appears  and,  once  the 
beep  sounds,  the  "Program  Master  Menu"  is  ready  for  use.  This  "page" 
is  shown  at  Fig.  3.  On  this  menu,  the  user  can  select  from  six 
options.  These  options  include: 

(1)  Quit  (exit  the  program), 

(2)  To  change/view  program  parameters, 

(3)  To  execute  the  graphics  routine, 

(4)  To  print  out  a  table  of  data, 

(5)  To  run  the  simulation, 

(6)  To  determine  pressure  losses. 
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Should  Che  user  select  the  "change/view  program  parameters" 
option,  that  menu  will  be  the  next  on  the  terminal  screen.  Fig.  4 
shows  an  example  of  this  menu.  On  this  menu  the  user  again  has 
several  options.  These  include: 

(1)  Quit  (exit  the  program), 

(2)  To  change  drilling  fluid  parameters, 

(3)  To  change  drill  pipe  parameters, 

(4)  To  change  wellbore/formation  parameters, 


(5)  To  change  casing  parameters, 


(6)  To  return  to  the  main  program  menu. 


Each  of  these  menus  will  then  take  the  user  to  the  appropriate  menu  to 


change  and/or  view  the  program  parameters.  If  the  user  wishes  to 


change/view  the  drilling  fluid  parameters,  he  would  press  the 


appropriate  key  and  the  menu  would  be  revealed  on  the  terminal  (Fig. 


5).  The  drilling  fluid  parameters  built  into  the  program  include 


thermal  conductivity  of  the  mud,  inlet  mud  temperature,  mud  density, 


flow  rate,  viscosity  of  the  mud,  and  specific  heat  of  the  mud.  Also, 


there  is  an  option  to  return  to  the  drilling  parameters  menu. 


If  the  user  had  desired  to  change/view  one  of  the  other  parameter 


menus,  a  similar  operation  could  be  conducted.  Each  of  these  other 


menus  would,  again,  lead  the  user  through  the  required  entries,  and 


back  to  the  main  program  menu.  Figs.  6  through  8  show  these  menus. 


Once  the  user  is  satisfied  with  the  values  of  each  of  the  parameters, 


other  operations  can  be  performed. 


TV 
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The  quickest  and  easiest  way  to  understand  data  is  by  visual 
inspection,  either  in  the  form  of  pictures  or  graphical  plots.  Both 
of  these  method  are  used  in  the  heat  transfer  model.  Finally,  the 
results  of  the  simulation  can  be  saved  on  a  disk  and/or  routed  to  a 
printer.  With  these  options,  the  presentation  and  understanding  of 
the  simulation  can  be  maximized. 

MODEL  FORMULATION 

The  mathematical  system  of  equations  used  in  this  simulator  are 
from  the  papers  by  Ramey®  and  Holmes  and  Swift^®.  Ramey’s  solution 
method  was  further  refined  and  amplified  by  Horne  and  Shinohara^®. 
Incropera  and  DeWitt^  and  Wilhite's  paper^®  were  used  as  the  basis 
for  the  calculation  of  the  over-all  heat  transfer  coefficients,  U. 

Generally,  there  are  three  components  of  wellbore  thermal 
simulators:  flowing  stream,  well  completion,  and  formation.  The 

relative  importance  of  each  component  depends  on  the  particular 
application  and  the  desired  result.  For  predicting  surface 
temperatures  in  high-rate  production  wells,  formation  calculation  is 
much  less  important  than  flowing  stream  calculation^.  This  is 
particularly  applicable  in  geothermal  wells.  On  the  other  hand,  for 
predicting  permafrost  thawing  around  an  arctic  well  the  flowing  stream 
is  less  important  than  the  formation. 

Circulation  of  fluid  during  the  drilling  operation  is  represented 
schematically  in  Fig.  9.  The  process  of  circulation  has  three 
distinct  phases^:  (1)  fluid  enters  the  drill  pipe  at  the  surface  and 

passes  down  the  drill  pipe;  (2)  fluid  exits  the  drill  pipe  through  the 
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Fie.  9  -  Schematic  of  circulating  fluid  system. 


bit  and  enters  the  annulus  at  the  bottom  of  the  wellbore;  and  (3) 
fluid  passes  up  the  annulus  and  exits  the  annulus  at  the  surface.  To 
simulate  the  thermal  behavior  in  the  system,  each  of  the  phases  of 
circulation  must  be  described  mathematically. 

In  Phase  1,  the  fluid  enters  the  drill  string  at  a  specified 
temperature,  TDo .  As  the  fluid  passes  down  the  pipe,  its  temperature 
is  determined  by  the  rate  of  heat  convection  down  the  drill  pipe,  the 
rate  of  heat  exchange  between  the  drill  pipe  and  the  annulus,  and 
time.  Phase  2  of  the  circulating  operation  simply  requires  that  the 
fluid  temperature  at  the  exit  of  the  drill  string  be  the  same  as  the 
fluid  temperature  at  the  entrance  of  the  annulus;  i.e.,  Tj)(Z,t)  - 

1  O 

T^(Z ,  t)  .  Therefore  in  Phase  3,  the  fluid  enters  the  annulus  at 
TD(Z,t).  As  the  fluid  flows  up  the  annulus,  its  temperature  is 
determined  by  the  rate  of  heat  convection  up  the  annulus,  the  rate  of 
heat  exchange  between  the  annulus  and  the  drill  pipe,  the  rate  of  heat 
exchange  between  the  formation  adjacent  to  the  annulus  and  the  fluid 
in  the  annulus,  and  time.  These  rates  of  heat  exchange  and  the  time 
dependency  of  the  mud  temperature  are  described  by  well  known  heat- 

*\  c 

flow  equations.  J  Consequently,  as  shown  in  Appendix  A,  the 
temperature  profiles  in  the  drill  pipe,  annulus,  and  formation  can  be 
obtained  by  solving  Eqs .  1  through  3  with  the  added  condition  as 
stipulated  in  Eq.  4,  once  appropriate  initial  and  boundary  conditions 
are  specified. 

Eq .  1  describes  the  heat  flow  within  the  drill  string.  The  terms 
on  the  left  hand  side  represent  the  vertical  and  radial  convective 


heat  transfer,  respectively.  The  terms  on  the  right  represent  the 
accumulation  of  energy  within  the  drill  string. 


ADpvDC  5Tp(Z,t)  +  2*rDU[TD(Z,t) 


TA(Z,t)]  -  -pADCp  5TP(Z|t) 


The  energy  balance  in  the  flowing  annulus  is  represented  by  Eq. 
2.  The  three  terms  on  the  left  hand  side  represent  the  vertical 
convective  heat  transfer  within  the  fluid,  radial  convection  between 
the  drilling  fluid  and  the  drill  pipe  wall,  and  radial  convection 
between  the  drilling  fluid  and  the  casing  or  formation,  respectively 
Heat  accumulation  is  accounted  for  by  the  term  on  the  right. 


AAPVAcp 


«TA(Z,t) 


+  2*rDU[TD(Z,t) 


-  TA(Z,t)]  +  Zwrghf  [Tf(rw,  Z,  t) 

-  TA(Z,t>]  -  pAnC  6TA(Z>t) 

St 


Eq.  3  is  a  two-dimensional  thermal  conductivity  equation 
representing  heat  flow  in  the  formation.  The  terms  on  the  right 
account  for  the  vertical  and  radial  conduction,  and  the  term  on  the 
left  accounts  for  the  heat  accumulation. 


5Tf(rw,Z,t)  kf  1  S  I"  STf(rw,Z,t) 


St  p fcpf  r  6 r 


'V'  «. 

■  v  -  --.  s.  S. 


r  6  r 


One  boundary  condition  for  Eq.  3  requires  that  the  flux  out  of 
the  formation  be  the  same  as  the  flux  into  the  annulus  and  is  given  by 
Eq .  4.  Since  at  large  values  of  r  the  geothermal  temperature  is 
undisturbed,  the  second  boundary  condition  is:  Tf(r-*»,Z,t)  -  Ts  +  GZ. 
Thus , 


2jrrBhf[Tf(Z,t)  -  TA(Z,t); 


2irrBkf 


S Tf(rw,Z,t) 


Sr 


(M 


rB 


To  obtain  the  energy  balances  above  which  describe  the  thermal 
behavior  of  the  wellbore  certain  assumptions  about  the  heat  transfer 
mechanisms  and  flow  behavior  are  required^: 

1.  Flow  is  steady  state  and  fully  developed. 

2.  Flow  is  turbulent  in  the  drill  pipe  and  across  the  drill  bit, 

and  laminar  in  the  annulus . 

I 

3.  Heat  transfer  within  the  drilling  fluid  is  by  axial 
convection.  Axial  conduction  is  neglected. 

4.  The  radial  temperature  gradient  within  the  drilling  fluid  may 
be  neglected. 

5.  Heat  generation  by  viscous  dissipation  within  the  fluid  may 
be  neglected. 

i 

|  6.  Fluid  properties  such  as  density,  thermal  conductivity,  and 


specific  heat  are  independent  of  temperature. 


HOLMES  AND  SWIFT^  MODEL 


The  equations  in  this  portion  of  the  simulator  are  based  on  the 
assumptions  that  the  heat  transfer  between  the  fluid  in  the  annulus 
and  the  formation  can  be  approximated  by  steady-state  linear  heat 
transfer.  The  work  of  Edwardson  et  al . ^  has  shown  that  the 
temperature  is  relatively  constant  at  any  point  sufficiently  removed 
from  the  drill  bit.  This  effect  shows  that  the  steady-state 
assumption  appears  to  be  a  close  enough  approximation  of  this 
phenomenon.  Other  simplifying  assumptions  are  that  the  heat  generated 
by  the  drill  bit  is  negligible  and  that  a  linear  geothermal  profile 
exists . 

The  development  of  the  model  is  depicted  in  Fig.  10.  A  slab  of 
thickness  dx  is  used,  assuming  heat  transfer  in  the  radial  direction 
and  no  significant  longitudinal  conduction.  The  heat  accumulation  of 
the  fluid  in  the  annulus  between  depth  x  and  x  +  dx  is  given  by 


Qax  '  Qa(x+dx)  ”  mcp[^Ax  '  ^A(x+dx) 1 


(5) 


and  the  steady-state  approximation  of  the  heat  transferred  between 
the  annular  fluid  and  the  formation  is  given  by 


*  -  2*rwU(TA  -  Tf)dx 


(6) 


The  heat  balance  across  the  drill  pipe  is  represented  by 


Fig.  10  -  Differential  fluid  element  used  to  derive  the  Holmes 
and  Swift^  model. 


Qap  “  2^rj)hp(Tp 


TA)dx 


Combining  these  equations  yields  the  over-all  heat  transfer  through 
the  annulus , 

n>cp  ^  +  2jrrDhp(Tp  -  TA)  -  2*rwU(TA  -  Tf)  . 

dx 

The  formation  temperature  may  be  approximated  as: 

Tf  -  Ts  +  GZ  . 

Substituting  Eq  9  into  the  over- all  heat  balance  for  the  annulus 
produces  the  heat  balance  for  the  element  across  the  annulus  fluid: 

dT 

™-p  *  2»rDhp(Tp  -  Ta)  -  2*rwU(TA  -  Ts  -  GZ)  .  ( 

dx 

■'*.  similar  development  for  the  fluid  in  the  drillstem  gives  the 

following  heat  balance: 


P  -  2irrDhp  (Tp  -  Ta) 


( 


These  are  the  equations  of  the  linear  heat  transfer  model  of  Holmes 
and  Swift^V  The  equations  are  then  integrated  into  their 


applicable  form.  These  are,  for  the  temperature  of  the  mud  in  the 
drillstem, 


Tp  -  K1eClZ  +  K2eC2Z  +  GZ  +  Ts  -  GA  .  . 
and  for  the  temperature  of  the  mud  in  the  annulus, 

Ta  -  KLC3eClZ  +  K2C4eC2Z  +  GZ  +  Ts 

where 

Cx  -  (B/2A) [ 1  +  (1  +  4/B)° • 5 ]  . 

C2  -  ( B/2A) [ 1  -  (1  +  4/B)0'5]  . 

C3  -  1  +  B/2 [ 1  +  (1  +  4/B)° • 5 ]  . 

C4  -  1  +  B/2 [ 1  -  (1  +  4/B)0-5]  . 

A  -  mcp/2^rDhp  . 

B  -  rwU/rDhp  . 


These  equations  when  applied  with  the  proper  boundary  conditions 


represent  the  analytical  solution  of  the  mud  temperature  profiles  for 
the  fluid  in  the  drillstem  and  annulus. 

Derivation  of  Circulating  Temperature  Equations 

Since  the  annular  and  drill  pipe  mud  temperatures  are  equal  at 
the  bottom  of  the  wellbore,  the  following  boundary  conditions  may  be 
applied  to  obtain  the  bottomhole  temperature^. 


Boundary  Condition  1  at  x  -  0;  Tq  -  TDo 


I 


Boundary  Condition  2  at  x  -  L;  THp  -  THa 

For  these  boundary  conditions  the  following  integration  constants  are 
obtained: 


K1  -  TDo  '  K2  -  Ts  +  GA . (20) 

K?  -  GA  -  t TPo  -  Ts  +  GA]eCIL  (1  -  C3)  . (21) 

eC2L  (1  -  C4)  -  eClL  (1  -  C3) 


These  constants  of  integration  are  applied  to  Eqs .  12  and  13  in 
I  order  to  calculate  the  temperature  at  any  point  in  the  well  during 


circulation . 


RAMEY0  MODEL 


The  classic  study  performed  by  Ramey  on  wellbore  heat 

transmission  derived  the  temperature  distribution  in  a  well  used  for 

1  ft 

injecting  hot  fluids.0  Ramey  later  expanded  on  this  to  give  the  rate 

of  heat  loss  from  the  well  to  the  formation. ®®  However,  by  assuming 

that  the  fluid  remains  at  its  inflow  temperature,  Ramey's  analysis 

effectively  gave  the  heat  loss  at  infinite  fluid  flow  rate  --  in  other 

words,  the  maximum  possible  heat  loss  rate.-'-®  A  paper  presented  by 

Horne  and  Shinohara'-®  reexamined  this  problem  for  finite  fluid  flow 

and  determines  the  heat  loss  rate  as  a  function  of  fluid  properties 

and  fluid  flow  rate.  Their  paper  presents  results  for  both  producing 

and  injection  wells,  which  is  valuable  when  analyzing  geothermal 

wells.  It  is  assumed  the  only  single-phase  fluids  are  flowing  in  the 

wellbore.  For  single-phase  flow,  the  formulation  permits  direct 

calculation  of  wellbore  heat  loss  with  various  production  and 

1  ft 

injection  conditions.  °  Derivation  of  Ramey's  wellbore  heat 
transmission  solution  is  found  in  Appendix  g. 31-34 

Ini ection 

As  in  Fig.  11,  consider  a  heat  balance  in  the  radial  direction  on 
a  section  of  a  well  with  height  dz ,  losing  heat  at  rate  dq  from  the 
casing  to  the  formation.  Then, 


dq 

dz 


2,,kEr°"  <TD  -  Tt) 
kf  +  rDUf(t) 


(22) 


For  long  periods  of  time,  f(t)  can  be  approximated  as 


f(t)  -  -In  r2  -  0.290  +  0(r2/4c«t) 

2jat 


(23) 


Ramev's  Eq .  23  is  adequate  for  use  for  periods  of  time  greater 
than  one  week.  For  time  periods  less  than  one  week,  Ramey  provided  a 
graphical  solution.  However,  in  order  to  further  refine  the  computer 
model  to  more  efficiently  use  Ramey's  graph  of  the  function  f(t) 
versus  log^Q  (at/r2^),  the  curve  representing  constant  temperature  at 
r  -  r2 .  cylindrical  source,  was  digitized  using  the  computer  program, 


DIGIT,  which  is  available  on  the  Petroleum  Engineering  Department's 
PRIME  computer.  Once  the  curve  was  digitized,  the  data  were  then 
analyzed  by  the  PRIME'S  least  squares  curve  fitting  program,  CURVIT. 
It  was  determined  that  the  data  provided  by  Ramey's  curve  fit  a  third 
order  power  curve  (log-log)  almost  exactly.  Fig-  12  is  a  plot  of  the 
third  order  power  curve  generated  by  the  computer  versus  the  actual 
data  provided  by  Ramey.  The  curve  in  question  is  modeled  by  the 
equat ion 


f(t)  -  0.0018302  (In  x)3  -  0.095016  (In  x>: 


+  0.49045  (In  x)  -  0.44056 


Performing  an  over-all  heat  balance  on  the  well  and  considering 


the  changing  temperature  of  the  fluid  as  it  flows  in  the  tubing,  Tp 
can  be  estimated  as 


GZ  +  Ts  -  GA  +  (TDo 


Ts  +  GA) e 


■Z/A 


(25) 


where  GZ  +  Ts  is  Tj,  the  formation  temperature  (assuming  linear 
geothermal  gradient),  Ts  is  the  surface  temperature,  and  Z  is  measured 
downward.  Tqq  is  the  injection  temperature.  A  is  a  group  of 
variables  which  Ramey  defined  as 


A(t)  -  WCP[kf  *  Uf(t)1 
2jrkfrDU 


(26) 


Integrating  Eq .  22  (with  respect  to  depth  Z)  and  substituting  TD 
”  ^Do ■  Ramey  obtained 


q<n  _  2*rPUkf 

kf  +  rDUf(t) 


(TDo 


TS)L 


(27) 


as  the  total  heat  loss  rate  from  a  well  of  total  depth  L.  However,  bv 
substituting  the  actual  value  of  Tq  from  Eq .  25, 


q 


(TDo  +  GA  -  Ts )  (1 


(28) 


or  alternatively, 


q  -  WCP (tDo  -  Td)  . (29) 

Comparing  [the  total  heat  loss  rate  for  a  well  with  infinite 
flow  rate  (Eq.  27)]  with  q  [the  total  heat  loss  for  a  well  with  flow 
rate  w  (Eq.  28)],  q  is  found  to  be  a  function  of  both  w  and  Cp , 
whereas  q^  is  independent  of  both  the  fluid  that  flows  in  the  wellbore 
and  the  rate  at  which  it  flows. 

Production 

With  a  well  producing  hot  fluids,  Eqs .  25,  27,  and  28  can  be  used 
by  replacing  Ts  with  TDq.  With  a  geothermal  well  producing  hot  fluid 
without  tubing,  the  temperature  in  the  well  as  a  function  of  height  y 
above  the  producing  depth  can  be  given  as 


TD  “  (TDo  *  Gy)  +  GA(1  -  e'y/A) . (30) 

Here,  the  temperature  of  the  earth  is  Tjj0  -  Gy  and,  since  there  is  no 
tubing,  U  is  infinite;  therefore, 

A  -  WCP£(t) . (31) 

2?rkf 


Then,  the  maximum  total  heat  flow  rate  from  a  well  of  depth  L  would  be 


whereas  the  actual  total  heat  flow  rate  would  be 


q  =  Gwc  [  L  +  A(e  L//A-  1)  ] 


Clearly,  for  large  flow  rates,  the  heat  loss  rate  approaches  that 
for  infinite  f low^-® .  if  the  exponential  function  in  Eq.  33  (or  Eq. 

28)  is  expanded  in  a  polynomial  series,  the  significance  of  other 
governing  parameters  can  be  evaluated^-®.  Substituting 


i  -  ±_  * 

A  2A2  6A- 


then , 


7rHkfi/  wcnGLJ 
q  =  r  -  P  + 


c-Dn  wcPGLn 


Comparing  this  with  Eq .  32  for  q^,,  the  difference  between  q  and  q^,  (1) 
is  always  negative,  q  <  q^;  (2)  increases  with  geothermal  gradient,  a; 
(3)  increases  rapidly  with  depth,  L;  (4)  decreases  with  increasing 


rv 


flow  rate,  w  (since  A  is  a  linear  function  of  w  and  Cp) ;  (5)  decreases 
with  increasing  specific  heat,  Cp ;  for  example,  a  steam  well  will  lose 
less  heat  than  a  hot  water  well;  and  (6)  increases  with  time,  t  [since 
f(t)  increases  with  t  and  A  decreases]-'-®. 

Thus,  the  traditionally  used  formulas  for  heat  loss  from  a  well 
carrying  hot  fluid  in  single-phase  flow  are  shown  to  overestimate  heat 
loss  severely,  particularly  in  deep  wells  with  moderate  -  to  -  low  flow 
rates  of  low  specific-heat  fluids.  This  would  be  important  in  steam 
injection,  where  the  quantity  of  steam  required  would  be 
overestimated,  or  when  flow  testing  a  geothermal  well,  where  down-hole 
enthalpies  and  steam  quality  are  calculated  by  considering  wellhead 
conditions  and  wellbore  heat  loss.  A  more  complete  discussion  of  the 
over-all  heat  transfer  coefficient,  U,  can  be  found  in  Appendix  D. 


PRESSURE  LOSS  CALCULATIONS 


The  basis  of  any  simulation  of  the  drilling  process  is  the 
drilling  fluid  or  mud.  The  drilling  fluid  serves  a  variety  of 
purposes  such  as  cooling  the  bit,  transporting  the  cuttings  to  the 
surface,  and  controlling  subsurfac  '  pressures.  Drilling  fluids  have 
progressed  from  little  more  than  clay  suspensions  to  highly  complex 
substances  both  Theologically  and  chemically.  This  problem  is  further 
aggravated  by  the  fact  that  the  rheology  and  chemical  make  up  can 
change  greatly,  even  during  the  course  of  drilling  a  single  well. 

Knowing  what  displacement  pressure  and  flow  rate  will  maintain  a 
drilling  fluid  or  cement  slurry  in  turbulent  or  plug  flow  in  the 
wellbore  annulus  is  essential  in  the  design  of  a  drilling  hydraulics 
program  or  primary  cement  job.  The  frictional  pressure  loss  term  in 
the  pressure  balance  equation  is  the  most  difficult  to  evaluate^5. 

The  pressure  balance  equation  in  question  is 

Pl  +  0.052p(D2  -  D]_)  -  8.074xl0’4p(v22  -  vx2) 

+  App  -  Apf  -  p2 .  (35) 

Fluids  in  plug  or  turbulent  flow  exert  a  uniform  displacement  force 
against  the  mud  in  the  wellbore  annulus^  In  laminar  flow,  cement 
has  a  parabolic  or  "bullet-shaped"  velocity  profile  across  the  area  of 
flow.  This  results  in  the  cement  "jetting"  through  the  drilling 
fluid^-  Incomplete  mud  removal  can  result  in  poor  cement  bonding, 
zone  communication,  and  ineffective  stimulation  treatments. 
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A  mathematical  or  empirical  fluid  model  describes  the  flow 
behavior  of  a  fluid  by  expressing  some  type  of  relationship  between 
shear  rate  and  shear  stress .  This  model  of  the  viscous  forces  present 
in  a  drilling  fluid  is  required  for  the  development  of  frictional 
pressure  loss  equations.  For  a  Newtonian  fluid,  the  ratio  of  the 
shear  stress  to  the  shear  rate  is  a  constant.  However,  for  non- 
Newtonian  fluids  the  relationship  between  shear  stress  and  shear  rate 
is  more  complicated.  Most  drilling  fluids  are  non-Newtonian  in 
nature.  As  of  the  date  of  this  paper,  a  generalized  relationship  for 
all  non-Newtonian  fluids  has  not  been  found^ .  Instead,  several 
models  have  been  proposed  to  describe  the  behavior  of  various  ideal 
non-Newtonian  fluids. 

For  this  model,  three  different  rheological  models  were  chosen 
because  of  their  general  use  in  the  drilling  industry.  The  three  most 
commonly  used  and  understood  models  are  the  Newtonian,  Bingham 
Plastic,  and  the  Ostwald-deWaele  Pseudoplastic  (Power  Law).  Appendix 
C  contains  a  complete  discussion  of  the  equations  used  for  pressure 
loss  calculations  for  each  of  the  three  most  commonly  accepted  fluid 
rheo logies . 

NEWTONIAN  FLUID  MODEL 

The  Newtonian  fluid  model  is  defined  by  the  relationship 


r  -  H  7 


(36) 


At  conditions  of  constant  temperature  and  pressure,  the  shear 
rate  and  the  shear  stress  are  directly  proportional.  The  constant  of 
proportionality,  n ,  is  called  the  absolute  viscosity.  Fig.  13  shows 
the  flow  curves  of  Newtonian  fluids.  Note  that  the  curves  are 
straight  lines  which  pass  through  the  origin,  and  that  the  slopes  of 
the  lines  are  n ,  the  absolute  viscosity.  As  shown  in  Fig.  13,  the 
thicker  the  consistency  of  the  fluid,  the  larger  is  the  magnitude  of 
the  slope. 

Water  and  several  pure  organic  liquids  are  Newtonian  fluids  in 
nature-*'.  Drilling  fluids  rarely  behave  as  Newtonian  fluids.  Because 
Newtonian  fluids  behave  in  a  relatively  simple  manner,  they  provide  an 
ideal  fluid  for  fluid  flow  experiments®®. 


Non-Newtonian  Fluid  Models 

The  ratio  of  the  shear  rate  to  the  shear  stress  of  a  non- 
Newtonian  fluid  is  not  constant.  This  is  true  for  most  drilling 
fluids.  The  two  most  widely  accepted  mathematical  models  for 
describing  non-Newtonian  drilling  fluids  are  called  the  Bingham 
Plastic  model  and  the  Ostwald  -  deWaele  power  law  model.  Hence  it  is 
generally  accepted  that  drilling  fluids  are  typified  by  one  of  these 
two  models. 


Bingham  Plastic  Model 

The  Bingham  Plastic  model  is  defined  by  the  relationship 


’00 


As  can  be  seen  in  Fig.  14,  the  Bingham  fluid  differs  from  a  Newtonian 
fluid  most  notably  by  the  presence  of  a  yield  stress,  most  commonly 

-I 

called  the  "yield  point"J/.  The  yield  stress  is  a  measure  of  the 
electrical  attractive  force  in  the  drilling  mud  under  flowing 
conditions.  The  Bingham  Plastic  model  assumes  that  a  fluid  behaves  as 
an  elastic  solid  up  to  the  yield  stress^.  Then  it  behaves  as  a 
Newtonian  fluid  above  this  point.  No  bulk  movement  of  the  fluid 
occurs  until  the  applied  stress  exceeds  the  yield  stress.  Once  the 
yield  stress  is  exceeded,  equal  increments  of  shear  stress  produce 
equal  increments  of  shear  rate.  Forces  less  than  the  yield  stress 
produce  a  deformation,  usually  ignored  because  it  is  so  small,  but 
cause  no  discernable  flow.  Bingham  Plastic  -  type  fluids  flowing 
through  a  tube  are  characterized  by  an  inner  plug  moving  in  a  fluid 
ring.  Attempts  to  fit  Bingham  Plastic  fluids  with  a  power  law  -  type 
model  result  in  extremely  low  values  for  the  Power  Law  exponent,  n. 

The  apparent  viscosity,  or  effective  viscosity,  defined  as  the 
ratio  of  the  shear  stress  to  the  shear  rate,  varies  with  shear  rate 
for  non-Newtonian  fluids.  The  apparent  viscosity  is  the  slope  of  a 
line  from  the  origin  to  some  particular  shear  rate.  The  slopes  of  the 
dashed  lines  in  Fig.  14  represent  apparent  viscosities  at  various 
shear  rates.  The  apparent  viscosity  decreases  with  increased  shear 
rate.  This  phenomenon  is  called  "shear  thinning."  As  shear  rates 
approach  infinity,  the  apparent  viscosity  reached  a  limit  called  the 


SHEAR  RATE  (  1/SEC  ) 


48 


"plastic  viscosity" ^ .  The  plastic  viscosity  is  the  slope  of  the 
Bingham  Plastic  flow  curve. 

This  model  has  been  used  extensively  in  the  oil  industry.  The 
model  is  easy  to  use,  and  it  represents  many  drilling  fluids 
reasonably  well.  In  fact,  the  commonly  used  Fann  V-G  meter  was 
specifically  designed  to  facilitate  the  use  of  the  Bingham  Plastic 
model  in  the  field.  However,  the  Bingham  model  usually  does  not 
represent  drilling  fluids  at  low  shear  rates.  More  sophisticated 
fluid  models  are  gaining  a  wider  acceptance  in  the  oil  field  as 
engineers  become  more  familiar  with  their  use. 

Ostwald-deWaele  Pseudoplastic  (Power  Law-)  Model 

The  pseudoplastic  flow  curve  is  characterized  by  a  non-linear 
relationship  between  shear  rate  and  shear  stress  as  shown  in  Fig. 

15^.  Velocity  profiles  for  pseudoplastic  materials  tend  to  be 
flatter  than  those  for  Newtonian  fluids,  but  normally  a  plug  does  not 
develop-^.  The  relationship  for  the  typical  pseudoplastic  fluid  is 
given  by 

r  -  k  yn  .  (38) 

Generally  speaking,  the  consistency  factor,  k,  describes  the 
thickness  of  the  fluid  and  is  somewhat  analogous  to  apparent 
viscosity.  As  k  increases  the  mud  becomes  thicker.  The  flow  behavior 
index,  n,  indicates  the  degree  of  non-Newtonian  behavior.  When  n 


equals  one,  the  Power  Law  equation  becomes  identical  to  the  Newtonian 


as  dilatant;  the  apparent  viscosity  increases  as  shear  rate 
increases.  If  n  is  greater  than  zero,  but  less  than  unity,  the  fluid 
is  classified  as  pseudoplastic.  Pseudoplastic  fluids  exhibit  "shear 
thinning,"  which  means  that  the  apparent  viscosity  decreases  as  the 
shear  rate  increases.  For  drilling  purposes,  shear  thinning  is  a  very 
desireable  property,  and  the  majority  of  drilling  fluids  are 
pseudoplastic . 

The  Power  Law  model  is  widely  used  in  the  oil  industry  and  has 
replaced  the  Bingham  model,  to  some  extent.  The  Power  Law  model  is 
frequently  more  convenient  than  the  Bingham  model.  It  is  particularly 
suitable  for  graphical  techniques  since  rotary  viscometer  readings 
versus  rpm  and  flow- pressure  loss  versus  flow  rate  can  be  plotted  as 
straight  lines  on  log- log  paper.  The  Power  Law  model  also  more 
accurately  demonstrates  the  behavior  of  a  drilling  fluid  at  low  shear 


rates.  However,  the  Power  Law  does  not  include  a  yield  stress  and 
therefore  can  give  poor  results  at  very  low  shear  rates^ . 

Fig.  16  compares  the  flow  curve  of  a  typical  drilling  fluid  to 
the  flow  curves  of  Newtonian,  Bingham  Plastic,  and  Power  Law  fluids. 
The  typical  drilling  fluid  exhibits  a  yield  stress  and  shear  thinning. 
At  high  rates  of  shear,  all  models  represent  the  typical  drilling 
fluid  reasonably  well.  Differences  between  the  fluid  models  are  most 
pronounced  at  low  shear  rates.  The  Bingham  Plastic  fluid  includes  a 
simple  yield  stress,  but  does  not  accurately  describe  the  fluid 
behavior  at  low  shear  rates.  The  Power  Law  more  accurately  describes 
the  behavior  at  low  shear  rates,  but  does  not  include  a  yield  stress. 
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The  typical  drilling  fluid  exhibits  behavior  in  between  the  Bingham 
model  and  the  Power  Law  model. 

Recently  attempts  have  been  made  to  resolve  the  controversy  by 
proposing  models  that  combine  the  obvious  shear  thinning  features  of 
drilling  fluids  with  a  yield  stress.  These  proposed  models  include 
the  Casson,  Robertson  -  Stiff,  and  Herschel  -  Bulkley  models. 
However,  for  the  purposes  of  this  study,  these  models  will  not  be 
discussed. 


VERIFICATION  OF  THE  MODEL 


The  heat  transfer  and  temperature  distribution  prediction  model 
has  been  compared  with  data  presented  by  two  different  authors  to 
ensure  that  this  simulator  in  fact  does  present  reasonable  results. 

The  data  used  for  the  comparison  is  from  papers  by  Ramey®  and  Holmes 

1  O 

and  Swift1  .  Table  1  lists  the  information  used  to  compare  the  data 
generated  by  this  model  with  those  of  the  two  previously  mentioned 
references.  These  data  are  also  the  default  data  which  are  built  in 
to  the  program  and  from  which  all  parameter  variations  begin.  These 
data  are  actual  data  obtained  from  a  15,000  foot  Gulf  Coast  well. 

For  the  case  of  the  Holmes  and  Swift^®  paper,  the  well  was 
assumed  to  have  no  drill  bit  size  change  and  no  casing  set.  Fig.  17 
illustrates  the  temperature  profiles  of  the  mud  in  the  annulus  and  the 
mud  in  the  drill  pipe.  It  appears  that  the  maximum  mud  temperature 
generally  occurs  in  the  annular  fluid  at  some  point  above  the  bottom 
of  the  hole.  This  phenomenon  was  also  presented  by  the  Tragesser  et 
al .*-0  model  and  agrees  with  distributions  obtained  from  measured  mud 
temperatures  during  mud  circulation1-’.  For  the  purposes  of  providing 
a  beginning  point  for  Figs.  17  through  21,  the  assumed  geothermal 
profile  is  included  on  all  these  figures.  This  is  the  solid  curve  on 
each  graph.  Also,  for  simplicity,  only  the  annulus  temperatures  have 
been  plotted  in  the  figure.  However,  the  drill  pipe  temperatures 
obtained  for  each  run  would  follow  the  same  trend  as  that  shown  in 


TABLE  1  --  PROPERTIES  USED 

IN 

THE  COMPUTER 

MODEL  AND  THEIR  DEFAULT 

VALUES 

Property 

Units 

Default  Value 

DRILLING  FLUID: 

Thermal  conductivity 

Btu/(hr*f t 

*deg  F) 

1.0 

Inlet  temperature 

deg  F 

75.0 

Mud  density 

lb/gal 

10.0 

Flow  rate 

gal/min 

210.0 

Mud  viscosity 

lb/(ft*hr) 

110.0 

Specific  heat 

Btu/(lb*deg  F) 

0.4 

DRILL  PIPE: 

Heat  transfer 

Btu/(hr*sq 

ft 

coefficient 

*deg  F) 

30.0 

Inner  diameter 

inches 

5.965 

Outer  diameter 

inches 

6.625 

Incremental 

length  (delz) 

feet 

500.0 

WELLBORE/FORMATION : 

Depth 

feet 

15,000.0 

Temperature  gradient 

deg  F/ft 

0.0127 

Surface  air  temp 

deg  F 

70.0 

Diameter 

inches 

8.375 

Heat  transfer 

Btu/(hr*sq 

ft 

coefficient 

*deg  F) 

1.0 

CASING: 

Thermal  conductivity 

Btu/(hr*f t 

*deg  F) 

25.0 

Inner  diameter 

inches 

10.0 

Outer  diameter 

inches 

10.75 

Weight 

lb/ft 

51.0 

Grade 

N-80 

<  -“/-"I' 

GEOTHERMAL  (GRADIENT  =  0.0127  DEG  F/FT  FOR  ALL  FIGURES) 
DRILL  PIPE 


DEPTH  IN  THOUSANDS  OF  FEET 


PARAMETRIC  SENSITIVITY  ANALYSIS 


The  primary  purpose  in  performing  such  an  analysis  is  to 
establish  to  what  extent  variations  in  a  particular  parameter  will 
effect  the  temperature  distribution  in  the  wellbore.  This  is  of 
particular  interest  where  the  parameter  in  question  cannot  be  directly 
evaluated.  In  such  a  case,  the  degree  to  which  variations  in  this 
parameter  would  affect  the  temperature  indicate  how  important  is  an 
accurate  estimate  for  the  parameter.  It  also  provides  a  qualitative 
estimate  of  the  error  intrinsic  in  the  simulator  caused  by  using 
assumed  values. 

Some  parameters  which  have  a  distinct  influence  on  the 
temperature  profile  on  the  fluid  flowing  on  a  circulating  well  can  be 
directly  and  accurately  measured.  These  include  circulating  time, 
depth,  wellbore  geometry,  and  drilling  fluid  characteristics.  Hence, 
such  a  detailed  analysis  is  of  somewhat  less  importance  with  regard  to 
these  parameters,  but  nonetheless  appropriate  in  view  of  the  fact  that 
it  can  be  used  to  evaluate  the  validity  of  assuming  some  of  them  to 
remain  constant. 

Any  computer  simulator  of  such  an  imprecise  and  complicated  an 
operation  as  drilling  a  well  must  be  a  simplification.  The  fewer  the 
simplifications  and  assumptions,  the  more  accurate  the  model  will  be. 
Unfortunately,  there  is  no  physical  means  for  determining  the  dynamic 
temperature  distribution  in  a  well,  and  thus  there  is  no  independent 
method  of  accessing  the  accuracy  of  the  model.  For  this  reason,  any 
variation  in  a  parameter  that  caused  a  maximum  deviation  of  less  than 


5°F  was  considered  insignificant.  It  should  be  noted  that  interaction 
between  parameters  was  disregarded  and  that  if  several  of  the 
parameters  were  varied  simultaneously,  it  is  likely  that  the  total 
effect  could  be  significant,  whereas  individually  the  variation  had  an 
insignificant  effect. 

Fig.  18  demonstrates  the  effect  of  varying  wellbore  and  drill 
pipe  diameters  in  the  wellbore  on  the  temperature  of  the  mud  in  the 
annulus  during  circulation.  Table  2  and  3  provide  the  critical 
information  that  was  used  in  the  generation  of  the  data  for  the  plot 
in  Fig.  18.  Fig.  18  shows  that  as  the  diameter  of  the  wellbore  is 
increased,  the  bottomhole  temperature  also  increases.  However,  as  you 
decrease  the  diameter  of  drill  pipe,  the  bottomhole  temperature  tends 
to  decrease.  It  does  appear  that  this  effect  is  reversed  somewhat  in 
the  upper  portion  of  the  wellbore.  In  this  simulation,  the 
temperatures  in  the  annulus  are  approximately  equal  at  a  depth  of 
about  4,000  feet.  By  increasing  the  diameter  of  the  wellbore  to  12 
inches,  the  bottomhole  temperature  increased  about  8°F.  The 
bottomhole  temperature  was  found  to  decrease  by  about  16°F,  after 
decreasing  the  drill  pipe  diameter  to  4  inches. 

Fig.  19  illustrates  the  effect  of  varying  mud  circulation  rates. 
Table  4  lists  the  critical  data  that  was  used  in  the  generation  of  the 
data  for  the  plot  in  Fig.  19.  Fig.  19  reveals  that  as  the  mud 
circulation  rate  increases,  the  temperature  distribution  in  the 
annulus  decreases.  On  the  other  hand,  as  you  decrease  the  mud 
circulation  rate,  the  temperature  in  the  wellbore  increases.  A  33.3 


TABLE  2 


WELL  DATA  SUMMARY  FOR  VARYING  WELLBORE  DIAMETER 


Wellbore  diameter 

- 

varies  (inches) 

Drill  pipe  O.D. 

- 

6.625  (inches) 

Mud  flow  rate 

= 

300  (bbl/hr) 

Heat  transfer 
coefficient  (pipe) 

-  30 

(Btu/(hr*f t*°  F) ) 

Heat  transfer  -  1 

coefficient  (annulus) 

(Btu/(hr*ft*°F) ) 

Specific  heat 

0 

.4  (Btu/ (lb*°F) ) 

Mud  density 

- 

10.0  (lb/gal) 

Geothermal  grad. 

- 

0.0127  ( ° F/f t) 

Inlet  temperature 

- 

75  ( °F) 

Hydraulic  radius 

- 

0.4375  (inches) 

«•  «*’ 


TABLE  3 


WELL  DAT a  SUMMARY  FOR  VARYING  DRILL  PIPE  OUTER  DIAMETER 


Wellbore  diameter 

- 

8.375  (inches) 

Drill  pipe  O.D. 

- 

varies  (inches) 

Mud  flow  rate 

- 

300  (bbl/hr) 

Heat  transfer 
coefficient  (pipe) 

-  30 

(Btu/(hr*f t*°F) ) 

Heat  transfer  -  1 

coefficient  (annulus) 

(Btu/(hr*ft*°F) ) 

Specific  heat 

0 

.4  (Btu/( lb*° F) ) 

Mud  density 

- 

10.0  (lb/gal) 

Geothermal  grad. 

- 

0.0127  (°F/ft) 

Inlet  temperature 

mm 

75  ( ° F) 

Hydraulic  radius 

- 

0.4375  (inches) 

y 

Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


TABLE  4 


WELL  DATA  SUMMARY  FOR  VARYING  MUD  FLOW  RATE 


Wellbore  diameter 

— 

8.375  ( inches) 

Drill  pipe  O.D. 

- 

6.625  (inches) 

Mud  flow  rate 

- 

varies  (bbl/hr) 

Heat  transfer 
coefficient  (pipe) 

-  30 

(Btu/(hr*ft*°F) ) 

Heat  transfer  -  1 

coefficient  (annulus) 

(Btu/(hr*ft*°F) ) 

Specific  heat 

0 

.4  (Btu/(lb*°F) ) 

Mud  density 

- 

10.0  (lb/gal) 

Geothermal  grad. 

- 

0.0127  (°F/ft) 

Inlet  temperature 

- 

75  (°F) 

Hydraulic  radius 

- 

0.4375  (inches) 

3  Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


percent  decrease  in  the  flow  rate  resulted  in  a  20°F  increase  in  the 


bottomhole  temperature.  Increasing  the  flow  rate  by  33.3  percent 
resulted  in  a  17°F  decrease  in  the  bottomhole  temperature.  Although 
the  flow  rate  has  a  significant  effect  on  the  annular  temperature 
distribution,  it  can  be  reliably  measured.  As  was  the  case  in  Fig. 

18,  it  appears  that  this  effect  is  reversed  slightly  above  the  4,000 
foot  level  in  the  wellbore. 

Fig.  20  shows  the  effect  of  varying  the  heat  transfer  coefficient 
of  the  drill  pipe  on  the  temperature  of  the  annulus  during 
circulation.  Table  5  lists  the  data  used  for  the  generation  of  this 
plot.  As  demonstrated  in  Fig.  20,  as  the  heat  transfer  coefficient  of 
the  pipe  is  increased,  the  temperature  in  the  annulus  begins  to 
increase  also.  As  the  heat  transfer  coefficient  on  the  pipe  is 
decreased,  the  temperature  in  the  annulus  also  decreases.  An  increase 
in  the  heat  transfer  coefficient  of  66.7  percent  resulted  in  an 
increase  of  14°F  in  the  wellbore,  and  a  decrease  of  66.7  percent 
resulted  in  a  temperature  reduction  of  38°F.  This  is  obviously  a 
critical  parameter  to  be  considered  in  any  computer  model.  As  in 
previous  examples,  this  effect  also  appears  to  be  reversed  in  the 
upper  portion  of  the  well. 

Fig,  21  illustrates  the  effect  of  varying  the  heat  transfer 
coefficient  of  the  annulus  on  the  temperature  of  the  mud  in  the 
annulus  during  circulation.  Table  6  lists  the  data  that  was  used  to 
generate  the  data  for  this  plot.  As  shown  in  Fig.  21,  as  the  heat 
transfer  coefficient  in  the  annulus  is  increased,  the  temperature  in 


TABLE  5 


WELL  DATA  SUMMARY  FOR  VARYING  HEAT  TRANSFER 


COEFFICIENT  (PIPE) 


Wellbore  diameter 

- 

8.375  ( inches) 

Drill  pipe  O.D. 

- 

6.625  ( inches ) 

Mud  flow  rate 

- 

300  (bbl/hr) 

Heat  transfer 
coefficient  (pipe) 

-  varies  (Btu/(hr*ft*°F) ) 

Heat  transfer 
coefficient  (annulus) 

1  (Btu/(hr*f t*°F) ) 

Specific  heat 

- 

0.4  (Btu/(lb*°F) ) 

Mud  density 

- 

10.0  (lb/gal) 

Geothermal  grad. 

- 

0.0127  CF/ft) 

Inlet  temperature 

- 

75  ( ° F) 

Hydraulic  radius 

- 

0.4375  (inches) 

He 


u  Refer  to  the  appropriate  figure  to  obtain  information  regarding 
range  of  the  parameter  being  varied. 


TABLE  6  --  WELL  DATA  SUMMARY  FOR  VARYING  HEAT  TRANSFER 
COEFFICIENT  (ANNULUS) 


Wellbore  diameter  - 

Drill  pipe  O.D. 

Mud  flow  rate 

Heat  transfer  - 
coefficient  (pipe) 


8.375  ( inches ) 
6.625  (  inches ) 
300  (bbl/hr ) 
30  (Btu/(hr*ft*°F) ) 


Heat  transfer  -  varies  (Btu/(hr*ft*°F) ) 5 
coefficient  (annulus) 


Specific  heat  - 

Mud  density  - 

Geothermal  grad. 
Inlet  temperature  - 
Hydraulic  radius  - 


0.4  ( Btu/ ( lb*° F) ) 
10.0  (lb/gal) 
0.0127  ( ° F/f t ) 
75  ( ° F) 
0.4375  (inches) 


5  Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


V .  /■  V 


68 


Che  annulus  is  found  Co  also  increase.  As  Che  heac  cransfer 
coefficienC  of  che  annulus  decreases,  so  does  che  ceraperacure  of  Che 
mud  in  che  circulaCing  wellbore.  Increasing  che  heac  Cransfer 
coefficienC  by  100  percenC  resulced  in  a  CemperaCure  increase  of  10°F 
in  che  annulus.  A  decrease  of  100  percenC  resulced  in  a  CemperaCure 
decrease  of  approximately  20°F.  This  parameCer  is  difficulc  co  obcain 
an  accurace  escimace,  and  one  finds  much  wriCCen  in  che  liCeraCure 
regarding  ics  decerminacion. 

Fig.  22  shows  che  effecCs  of  varying  che  specific  heac  of  che  mud 
on  che  CemperaCure  of  che  annulus.  Increasing  or  decreasing  chis 
parameCer  by  50  percenC,  covers  Che  range  of  values  quoCed  in  che 
liCeraCure,  resulced  in  large  variacions  in  che  CemperaCure  of  che 
wellbore  annulus.  An  increase  of  50  percenC  in  Chis  parameCer 
resulced  in  a  boccomhole  CemperaCure  decrease  of  25°F.  A  similar 
decrease  in  Che  specific  heac  resulced  in  a  boccomhole  CemperaCure 
increase  of  30°F.  Since  a  small  change  in  che  value  of  chis  parameCer 
resulcs  in  a  large  alCeraCion  in  che  annular  CemperaCure  profile,  and 
as  che  variaCion  of  chis  parameCer  wichin  Che  syscem  cannoC  accuraCely 
be  assessed,  ic  is  very  significanC.  Table  7  liscs  a  summary  of  che 
combinacion  of  variables  used  for  chis  comparison. 

Increasing  che  fluid  densicy  from  10  lb/gal  Co  13.4  lb/gal 
decreased  che  boccomhole  CemperaCure  by  17°F,  while  increasing  ic  Co 
16  lb/gal  decreased  che  boccomhole  CemperaCure  by  29°F.  The  resulcing 
annular  CemperaCure  profiles  are  shown  in  Fig.  23.  Alchough  Che 
density  of  che  fluid  is  known  with  reasonable  accuracy  upon  entering 


TABLE  7  --  WELL  DATA  SUMMARY  FOR  VARYING  SPECIFIC  HEAT 


OF  THE  DRILLING  MUD 


Wellbore  diameter  — 

Drill  pipe  O.D. 

Mud  flow  rate  - 

Heat  transfer 
coefficient  (pipe) 

Heat  transfer  - 
coefficient  (annulus) 

Specific  heat  - 

Mud  density  =■ 

Geothermal  grad.  - 

Inlet  temperature  - 

Hydraulic  radius  - 


8.375  (inches) 
6.625  (inches) 
300  (bbl/hr) 
300  (Btu/(hr*ft*°F) ) 


1  (Btu/(hr*ft*°F) ) 


varies  (Btu/(lb*°F) )^ 
10.0  (lb/gal) 
0.0127  (°F/ft) 
75  ( ° F) 
0.4375  (inches) 


Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


Che  system,  its  variations  within  the  system  are  not.  In  particular, 
the  fluid  density  will  increase  due  to  the  addition  of  drill  cuttings 
during  drilling  operations.  Table  8  summarizes  the  variables  used  in 
this  particular  analysis. 

Alterations  in  the  geothermal  gradient  have  a  distinct  effect  on 
the  temperature  profile  in  the  annulus  as  illustrated  in  Fig.  24. 
Inasmuch  as  this  parameter  will  remain  constant  for  a  particular  well, 
it  cannot  be  determined  with  any  degree  of  precision.  Fig.  24 
indicates  that  a  poor  assumption  of  the  value  of  the  geothermal 
gradient  will  cause  errors  in  the  annulus  temperature  distribution. 

The  range  of  values  in  Fig.  24  is  representative  of  the  values  quoted 
in  the  literature.  Table  9  provides  a  list  of  the  parameters  used  in 
this  portion  of  the  study. 

The  inlet  fluid  temperature  can  be  accurately  measured.  However, 
Fig.  25  illustrates  the  need  to  insure  that  this  parameter  should  be 
continuously  monitored  during  drilling  operations  to  be  certain  that 
the  annulus  temperature  profiles  as  predicted  by  the  model  are 
accurate.  Variables  used  in  this  analysis  are  found  in  Table  10. 

Fig.  26  illustrates  that  the  temperature  profile  in  the  annulus 
is  relatively  insensitive  to  this  parameter  with  a  maximum  variation 
in  the  annulus  temperature  profile  of  about  20°F  at  the  largest 
hydraulic  radius  evaluated.  Table  11  list  the  values  of  the 
parameters  used  in  this  evaluation. 


These  data  and  graphs  agree  with  the  results  of  Holmes  and 
Swift^  very  closely.  When  comparing  Figs.  17  through  21  with  the 
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TABLE  8  --  WELL  DATA  SUMMARY  FOR  VARYING  MUD  DENSITY 


Wellbore  diameter  - 

Drill  pipe  O.D.  - 

Mud  flow  rate  — 

Heat  transfer  - 
coefficient  (pipe) 

Heat  transfer  - 
coefficient  (annulus) 

Specific  heat  - 

Mud  density  - 

Geothermal  grad.  - 

Inlet  temperature  - 

Hydraulic  radius  - 


8.375  (inches) 
6.625  (inches) 
300  (bbl/hr) 
30  (Btu/(hr*ft*°F) ) 


1  (Btu/(hr*f t*°F) ) 


0.4  (Btu/(lb*°F) ) 
varies  (lb/gal) ^ 
0.0127  ( °F/f t) 
75  ( ° F) 
0.4375  (inches) 


^  Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


TABLE  9 


WELL  DATA  SUMMARY  FOR  VARYING  GEOTHERMAL  GRADIENT 


Wellbore  diameter  - 

Drill  pipe  O.D.  - 

Mud  flow  rate 

Heat  transfer 
coefficient  (pipe) 

Heat  transfer  - 
coefficient  (annulus) 

Specific  heat  - 

Mud  density  - 

Geothermal  grad.  - 

Inlet  temperature  - 

Hydraulic  radius  - 


8.375  (inches) 
6.625  (inches) 
300  (bbl/hr) 
30  (Btu/(hr*ft*°F) ) 

1  (Btu/(hr*ft*°F) ) 

0.4  (Btu/(lb*°F) ) 
10.0  (lb/gal) 
varies  (°F/ft)8 
75  (°F) 
0.4375  (inches) 


O 

Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


inlet  fluid  temperature  on  the 


TABLE  10  --  WELL  DATA  SUMMARY  FOR  VARYING  INLET 


FLUID  TEMPERATURE 


Wellbore  diameter  - 

Drill  pipe  O.D.  - 

Mud  flow  rate  - 

Heat  transfer 
coefficient  (pipe) 

Heat  transfer 
coefficient  (annulus) 

Specific  heat  - 

Mud  density  - 

Geothermal  grad.  - 

Inlet  temperature  - 

Hydraulic  radius  - 


8.375  (inches) 
6.625  (inches) 
300  (bbl/hr ) 
30  (Btu/(hr*f t*°F) ) 


1  (Btu/(hr*ft*°F) ) 


0.4  (Btu/(lb*°F) ) 
10.0  (lb/gal) 
0.0127  ( °F/f t) 
varies  (°F)9 
0.4375  (inches) 


9  Refer  to  the  appropriate  figure  to  obtain  information  regarding 
the  range  of  the  parameter  being  varied. 


DEPTH  IN  THOUSANDS  OF  FEET 
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TABLE  11  --  WELL  DATA  SUMMARY  FOR  VARYING  HYDRAULIC  RADIUS 


■\ 


I 


Wellbore  diameter  - 

Drill  pipe  O.D. 

Mud  flow  rate  - 

Heat  transfer  - 
coefficient  (pipe) 

Heat  transfer  - 
coefficient  (annulus) 

Specific  heat  - 

Mud  density  - 

Geothermal  grad. 

Inlet  temperature  - 

Hydraulic  radius  - 


8.375  ( inches ) 
6.625  (inches) 
300  (bbl/hr) 

30  (Btu/(hr*ft*°F) ) 

1  (Btu/(hr*ft*°F) ) 

0.4  (Btu/( lb*° F) ) 
10.0  (lb/gal) 
0.0127  (°F/ft) 

75  (°F) 

varies  (inches)10 


10  Refer  to  the  appropriate  figure  to  obtain  information 
regarding  the  range  of  the  parameter  being  varied. 
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figures  of  Holmes  and  Swift  (Figs.  2  through  6)1  ,  very  similar 
results  can  be  obtained.  Therefore,  Figs.  22  through  26  are  most 
likely  to  be  representative  of  temperature  effects  on  the  annulus. 

Thus  it  has  been  shown  that  a  steady-state  analytical  model  can  be 
developed  to  determine  circulating  mud  temperature  distributions  for 
the  annulus  and  drill  string  fluids.  Figs.  27  through  30  provide  a 
comparison  of  the  data  provided  by  the  simulator  described  in  this 
thesis,  and  the  data  provided  by  holmes  and  Swift^®.  As  can  be  seen, 
excellent  agreement  exists.  Therefore,  this  model  provides  a  rapid 
method  for  computing  mud  temperature  profiles  and  bottomhole  mud 
temperatures  so  that  appropriate  mud  an  cement  properties  can  be 
selected. 

As  for  verification  of  the  Ramey®  model,  it  was  found  that  the 
computer  program  generated  very  good  agreement  with  Ramey's  results. 
However,  direct  comparison  with  any  data  generated  by  Ramey  could  not 

O 

be  accomplished.  Ramey0  presented  a  method  for  the  solution  of  the 
wellbore  heat  transfer  problem  which  was  a  combination  of  hand 
calculations  and  application  of  his  time  function  graph.  He  provides 
one  calculation  of  temperature  in  a  wellbore  at  a  certain  depth.  When 
this  point  was  compared  with  the  data  generated  by  the  computer  model, 
good  agreement  was  obtained.  However,  it  should  be  noted  that  Ramey's 
calculations  involved  round  off  error  which  was  not  present  in  the 
computer  model.  Fig.  31  is  a  graph  of  the  computer-generated  Ramey 
data.  It  also  provides  an  example  of  the  type  of  plotting  function 
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results  of  this  simulator  and  those 
same  mud  flow  rate.  200  BBL/HR. 
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which  is  available  in  the  computer  model.  This  plot  was  produced  by 
using  the  default  parameters  as  found  in  Table  1. 

It  is  also  possible  to  calculate  how  the  mud  temperature  profile 
is  affected  by  the  change  in  wellbore  sizes  due  to  drill  bit  changes 
during  various  phases  of  a  drilling  schedule.  Using  a  14,000  foot 
well  as  an  example,  it  will  be  assumed  that  there  was  one  change  in 
drill  bit  size  at  11,500  feet  and  changes  in  mud  circulation  rate  at 
10,000  feet  and  at  12,500  feet.  To  observe  these  effects,  the  mud 
temperature  profiles  would  be  calculated  for  the  mud  in  the  annulus 
and  in  the  drill  pipe  for  each  segment  of  drilling.  Fig.  32 
illustrates  a  part  of  each  mud  circulation  profile  for  each  additional 
segment  drilled.  The  end  of  each  segment  represents  the  bottomhole 
temperature  at  that  point  in  the  drilling  schedule.  By  constructing  a 
line  through  the  bottomhole  temperature  calculated  for  each  segment, 
we  can  observe  any  abnormalities  originated  by  each  change  in  the 
drilling  of  the  well. 

The  model  also  works  well  in  the  calculation  of  pressure  losses 
in  the  system  for  each  of  the  three  mud  rheologies  discussed: 
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Newtonian;  Bingham  plastic,  and;  power  law.  The  numbers  provided  by 
the  computer  model  were  verified  by  hand  calculations  and  by 
comparison  with  example  provided  oy  Bourgoyne  et  al . ^ . 

Representative  computer  output  for  the  three  mud  rheologies  is 
provided  in  Figs.  33  through  35. 

There  is  no  evidence  from  this  research  to  indicate  that  any  of 
the  major  assumption  enumerated  in  the  Introduction  are  invalid. 


Given  the  physical  geometry  of  the  system  being  modeled,  they  are  all 
logical  assumptions  and  the  problem  is  greatly  simplified  as  a  result. 

However,  the  last  assumption,  namely  that  the  fluid  properties 
are  independent  of  temperature,  is  open  to  debate.  It  is  evident  from 
Fig.  22  and  23  that  the  mud  specific  heat  and  mud  density  have  a 
significant  effect  on  the  temperature  distribution.  No  investigation 
was  conducted  to  assess  the  extent  to  which  these  parameters  are 
affected  by  temperature  and  until  this  is  done,  the  validity  of  this 
assumption  must  remain  in  doubt. 

APPLICATION  TO  FIELD  PROBLEMS 

The  mud  circulation  temperatures  obtained  by  this  method  have 
been  used  to  predict  logged  bottomhole  mud  temperatures.  Temperatures 
from  the  circulation  period  are  used  as  the  initial  temperatures  for 
inclusion  into  a  transient  line-source  solution  to  predict  the 
transient  mud  temperature  buildup  curves^.  Since  the  heat  flux  to 
the  formation  is  not  known  with  any  accuracy,  this  solution  is  an 
approximation  of  the  actual  temperature  buildup.  If  we  calculate  an 
approximate  heat  flux  by  using  Eq.  6,  the  solution  of  the  line-source 
equation  produces  a  mud  temperature  at  any  time  mud  circulation  is 
stopped.  In  most  cases,  where  the  diameter  of  the  wellbore  is  less 
than  a  foot,  the  mud  temperature  is  predicted  to  return,  within 
several  hours  after  mud  circulation  has  ceased,  to  approximately  the 
temperature  of  the  formation  in  the  vicinity  of  the  wellbore^. 

Data  for  the  15,000  foot  well  used  in  this  section  to  calculate 
the  temperature  profiles  as  shown  in  Fig.  17  were  applied  with  actual 


drill  bit  size  changes  made  during  the  drilling  of  this  well.  The 
bottomhole  temperatures  for  the  various  segments  drilled  were 
calculated.  These  bottomhole  temperatures  were  then  used  in 
conjunction  with  the  line -source  solution  for  the  temperature  buildup 
after  circulation  is  stopped.  The  predicted  logged  bottomhole  mud 
temperatures  to  be  obtained  at  the  time  of  logging  are  given  in  Table 
12.  The  actual  bottomhole  mud  temperature  are  also  presented.  Fig. 
36  shows  good  agreement  between  the  calculated  and  the  actual 
bottomhole  mud  temperatures.  Calculations  for  other  wells  have  shown 
similar  results,  which  shows  the  capability  of  the  model  to  simulate 
the  apparent  physical  phenomena^-* . 


TABLE  12  -  CALCULATED  AND  LOGGED  BOTTOMHOLE 
TEMPERATURES  FOR  A  15,000  FOOT  GULF  COAST  WELL 


Depth 

(ft) 

Wellbore  size 
(ft) 

Logged 

Bottomhole 

Mud  Temp . 

(°F) 

Calculated 
Bottomhole 
Mud  Temp . 
(°F) 

7,000 

0.771 

163 

160 

10,000 

0.508 

180 

187 

CONCLUSIONS 


The  primary  objective  of  this  work  was  to  develop  a  computer 
model  that  estimates  the  temperature  profile  in  a  wellbore  both 
accurately  and  quickly.  This  model  was  then  used  to  investigate  the 
significance  of  the  parameters  in  terms  of  their  effect  on  temperature 
distribution.  The  major  conclusions  of  this  work  are  as  follows: 

1.  The  use  of  interactive  techniques  and  graphics  greatly 
enhances  the  effectiveness  and  user  understanding  of  computer 
simulation. 

2.  Although  heat  transfer  and  temperature  distribution 
simulation  has  traditionally  been  left  to  those  with  access  to  large, 
high-speed  digital  computers,  the  personal  computer  is  well  suited  for 
applications  in  this  area. 

3.  It  has  been  shown  that  a  computer  simulation  using  a  steady- 
state  analytical  model  can  be  developed  to  determine  circulating  mud 
temperature  profiles  for  the  fluids  in  the  drill  pipe  and  in  the 
annulus . 

4.  This  simulator  can  provide  a  rapid  method  for  computing  mud 
temperature  profiles  and  bottomhole  mud  temperatures  so  that 
appropriate  drilling  mud  and  cement  properties  can  be  selected. 

5.  The  simulator  has  been  compared  with  actual  well  logging  data 
and  has  successfully  predicted  logged  bottomhole  temperatures. 

6.  The  following  parameters  have  a  crucial  effect  on  the 
wellbore  temperature  distribution: 

a.  drilling  fluid  flow  rate, 
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b.  heat  transfer  coefficient  of  the  drill  pipe, 

c.  heat  transfer  coefficient  of  the  annulus, 

d.  drilling  fluid  specific  heat, 

e.  drilling  fluid  density, 

f.  geothermal  gradient, 

g.  drilling  fluid  inlet  temperature, 

h.  depth. 

7.  The  simulator  has  demonstrated  that  large  errors  in  wellbore 
temperature  estimates  can  be  made  if  drilling  personnel  assume  that 
the  temperature  distribution  in  the  annulus  is  linear  in  nature. 

These  errors  can  range  as  high  as  nearly  40°F  at  the  boom  of  the  well, 
and  could  have  disastrous  effects  on  a  drilling  mud  program  or  cement 
job . 
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NOMENCLATURE 

2 

Aa  -  cross-sectional  area  of  annulus,  in^ 

•  2 

Aj)  -  cross-sectional  area  of  drill  pipe,  m4 
cl>c2>c3-c4  ”  Holmes  and  Swift13  equation  coefficients 
Cp  -  heat  capacity  of  fluid,  Btu/(lb-°F) 

Cp£  =  heat  capacity  of  formation,  Btu/(lb-°F) 

D  -  depth,  ft 
E  -  internal  energy 

f(t)  -  Ramey's®  transient  heat  conduction  time 
function  for  earth,  dimensionless 
g  -  gravitational  acceleration,  32.2  ft/sec 
G  -  geothermal  gradient,  °F/ft 
gc  -  conversion  factor,  32.2  ft- lb  mass/sec2-lb  force 
hf  -  wellbore  wall  heat  transfer  coefficient, 

Btu/(f t2- °F-hr) 

hp  -  over-all  heat  transfer  coefficient  across 
drill  pipe,  Btu/(ft2- °F-hr) 

H  -  enthalpy,  Btu/lb  mass 

J  —  mechanical  equivalent  of  heat,  778  ft-lb/Btu 
K  -  consistency  index,  eq.  cp 
,K2  -  Holmes  and  Swift13  integration  constants 
kf  -  formation  thermal  conductivity, 

Btu/(ft- °F-hr) 

L  -  total  depth  of  well,  ft 
ra  -  mass  flow  rate,  lb/hr 


wv-VAy.VtV. 


flow  behavior  index,  dimensionless 
pressure,  lb/in 

Ramey's®  heat  transfer  rate,  Btu/day 

heat  flow  in  the  annulus,  Btu/hr 

heat  flow  across  drill  pipe,  Btu/hr 

heat  flow  in  the  drill  pipe,  Btu/hr 

radial  space  variable,  in 

borehole  radius,  in 

drill  pipe  radius,  in 

outside  radius  of  casing,  ft 

wellbore  radius,  in 

annular  temperature,  °F 

drill  pipe  temperature,  °F 

formation  temperature,  °F 

inlet  drill  pipe  temperature,  #F 

bottomhole  temperature  of  fluid  in  annulus, 

°F 

bottomhole  temperature  of  fluid  in  drill 
pipe,  °F 

drill  pipe  temperature,  °F 
temperature  of  earth's  surface,  °F 
fluid  velocity 

over-all  heat  transfer  coefficient  between 
drill  pipe  and  annulus,  BTU/(ft^- °F-hr) 
average  flow  rate,  ft/sec 
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KV- 


V  -  specific  volume 
vA  -  annular  fluid  velocity,  ft/sec 
vD  -  drill  pipe  fluid  velocity,  ft/sec 
w  -  fluid  flow  rate ,  lb/day 
Wf  -  flow  work,  ft-lb  force,  lb  mass 
Z  -  depth  variable,  ft 
t  -  time  variable,  hr 

a  -  thermal  diffusivity  of  earth,  sq  ft  /  day 
(a  -  k/pcf) 

Ap  -  pressure  loss,  lb/in^ 

At  -  time  step  variable,  hr 

"  Tbottomhole  fluid  '  Toutlet>  °F 
p  -  fluid  density,  lb/gal 

P f  -  formation  density,  lb/ft-* 
p  -  fluid  viscosity,  cp 
Pco  -  plastic  viscosity,  dyne-sec/cm^ 

$  -  heat  flux  to  formation,  Btu/hr 
r  -  shear  stress,  dyne/cm^ 
rQ  -  yield  point,  lb/100  ft^ 

7  -  shear  rate,  sec"^ 
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APPENDIX  A 


MATHEMATICAL  MODEL  AND  TREATMENT 


Discussion  of  Mathematical  Model 

The  circulation  of  fluid  during  the  drilling  operation  has  been 
represented  schematically  in  Fig.  1,  and  the  process  of  circulation 
has  been  separated  into  three  well-defined  phases^ ;  (l)  fluid  enters 

the  drill  pipe  at  the  surface  and  passes  down  the  drill  string;  (2) 
fluid  exits  the  drill  string  through  the  nozzles  in  the  drill  bit  and 
enters  the  annulus  at  the  bottom;  and  (3)  fluid  passes  up  the  annulus 
and  exits  at  the  surface. 

The  following  assumptions  were  made  by  Raymond^  to  obtain  the 
desired  differential  equations: 

1.  Axial  conduction  of  heat  in  the  fluid  is  insignificant 
compared  with  axial  convection.  This  appears  to  be  an  excellent 
assumption  considering  what  normal  circulation  rates  are  experienced 
in  the  circulation  of  drilling  fluids. 

2.  There  are  no  radial  gradients  in  the  fluids  in  either 
the  drill  string  or  the  annulus.  Since  it  is  assumed  that  the  fluid 
in  the  drill  pipe  is  in  turbulent  flow  and  well  mixed  every  30  feet  in 
the  annulus  because  of  the  presences  of  tool  joints  and  drill  collars, 
this  is  also  an  excellent  assumption. 

3.  The  properties  of  the  fluid,  such  as  heat  capacity, 
density,  and  thermal  conductivity,  do  not  change  significantly  with 
increases  in  temperature. 


4.  Heat  generation  by  viscous  dissipation  in  the  fluid  is 


negligible . 

In  the  first  phase,  the  fluid  enters  the  drill  pipe  at  a 
specified  temperature,  Tj)0 .  As  the  fluid  passes  down  the  pipe,  its 
temperature  is  determined  by  the  rates  at  which  (1)  heat  is  convected 
down  the  drill  pipe,  (2)  heat  is  exchanged  between  the  drill  pipe  and 
the  annulus,  and  (3)  the  temperature  of  the  drill-pipe  fluid  changes 
with  time.  Consequently,  the  equation  that  describes  the  temperature 
of  the  fluid  in  the  drill  pipe  as  a  function  of  time,  t,  and  depth,  Z, 
is 


ADPvDC  6Tp(Z,t)  +  2*rDU[TD(Z,t) 
6  Z 


-  TA(Z,t>]  -  -pADC  *Tp(Z,t:) 

St 


(A-  1) 


In  Eq.  A-l,  it  has  been  assumed  that  the  over-all  heat  transfer 
coefficient,  U,  is  independent  of  depth  and  time.  Since  it  has 
already  been  assumed  that  the  fluid  thermal  properties  are  independent 
of  temperature,  this  is  a  good  assumption.  The  second  phase  of  the 
circulating  process  simply  requires  that  the  fluid  temperature  at  the 
exit  of  the  drill  pipe  be  the  same  as  the  temperature  at  the  entrance 
of  the  annulus;  i.e.,  TD(L,t)  -  TA(L,t).  Accordingly,  in  the  third 
phase,  the  drilling  fluid  enters  the  annulus  at  a  temperature  of 
Tjj(L,t).  As  the  fluid  is  pumped  up  the  annulus,  it  temperature  is 
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determined  by  the  rates  at  which  (1)  heat  is  convected  up  the  annulus, 
(2)  heat  is  exchanged  between  the  annulus  and  the  drill  pipe,  (3)  heat 
is  exchanged  between  the  formation  adjacent  to  the  annulus  and  the 
fluid  in  the  annulus,  and  (4)  the  temperature  of  the  annular  fluid 
changes  with  time.  Consequently,  the  equation  that  describes  the 
temperature  of  the  fluid  in  the  annulus  as  a  function  of  t  and  Z  is 

Aa PVAC  *TA(Z>t)  +  2*rDU[TD(Z,t) 

62 

-  Ta(Z, t) ]  +  2*rBhf[Tf(rw,Z,t) 

-  TA(Z , t) ]  -  pAnCpTA(Z’t)  .  (A-2) 

6 1 

Since  the  thermal  conductivity  of  the  formation  adjacent  to  the 
wellbore  is  small,  it  will  be  assumed  that  there  is  no  transfer  of 
heat  by  conduction  in  the  vertical  direction  in  the  formation.  Hence, 
the  formation  temperature  Tf(rw,Z,t)  is  regulated  by  the  radial 
diffusivity  equation 

6Tf(rw,Z,t)  kf  16 
5t  ^f^pf  £  6r 

Eq .  A- 3  is  coupled  with  Eq.  A-2  through  the  rate  of  heat  transfer 
between  the  fluid  in  the  annulus  and  the  formation.  This  requires 
that  the  flux  of  heat  out  of  the  formation  be  the  same  as  the  flux  of 
the  heat  into  the  annulus;  that  is,  at  any  depth  Z, 


r  5Tf(rw,Z,t) 


6r 


(A-3) 
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2*rBhf[Tf(Z,t)  -  TA(Z,t) 
27rrBkf 


STf(rw,Z,t) 


Thus,  the  temperature  profiles  in  the  drill  pipe,  annulus,  and 
formation  can  be  obtained  by  solving  Eqs.  A-l  through  A-3,  with  the 
additional  requirement  specified  by  Eq.  A-4  once  appropriate  initial 
and  boundary  conditions  are  specified. 


APPENDIX  B 


DERIVATION  OF  RAMEY'S8 
WELLBORE  HEAT  TRANSMISSION  SOLUTION 
In  1962,  Ramey  considered  the  injection  of  a  fluid  down  the 
tubing  in  a  well  which  is  cased  to  the  top  of  the  injection  interval. 
Assuming  that  the  fluid  is  injected  at  known  injection  rates  and 
temperatures,  he  determined  the  temperature  of  the  injected  fluid  as  a 
function  of  depth  and  time.  Fig.  B-l  presents  a  schematic 
illustration  of  the  problem.  As  shown  on  Fig.  B-l,  W  lb/day  of  fluid 
is  injected  in  the  tubing  at  the  surface  at  a  temperature  of  TQ.  The 
inside  radius  of  the  tubing  is  rj_,  and  the  temperature  T^  of  the  fluid 
in  the  tubing  is  a  function  of  both  depth  Z  and  t.  The  outside  radius 
of  the  casing  is  r^,  and  the  temperature  of  the  casing  outer  surface 
is  T2 ,  also  a  function  of  depth  and  time. 

The  usual  solution  procedure  for  flow  problems  of  this  type  is  to 
solve  the  total  energy  and  mechanical  energy  equations  simultaneously 
to  produce  both  temperature  and  pressure  distributions.  However,  the 
solution  may  be  approximated  by  the  following  considerations.  The 
total  energy  equation  is 

dH  -  g  dZ  +  ^  -  dQ  -  . (B-l) 

gcJ  ScJ  J 

Assuming  steady  flow  of  a  single  phase  fluid  in  a  pipe,  flow-work  Wf 
is  zero  and  Eq .  B-l  becomes 


If  it  is  assumed  that  the  fluid  flowing  is  an  incompressible 


i 

liquid,  the  kinetic  energy  term  becomes  zero.  Thus, 

dH  +  g  dZ  -  dQ . 

Sc  J 


(B-3) 


But  by  definition,  enthalpy  is 


dH  -  dE  +  d(PV) 
J 


dE 


+  V  dp 


(B-4) 


for  an  incompressible  liquid.  Or 


dH 


cdT 


V  dP 


(B-5) 


Neglecting  the  flowing  friction,  the  V  dP  term  is  equal  to  the  change 
in  fluid  head,  and  the  change  in  enthalpy  is 

dH  ~  cdT  +  g  dZ . (B-6) 

J 

Considering  flow  down  the  well,  the  increase  in  enthalpy  due  to 
increase  in  pressure  is  approximately  equal  to  the  loss  in  potential 
energy.  Conversely,  for  flow  up  the  well,  the  loss  in  enthalpy  due  to 
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che  decrease  In  pressure  is  approximately  equal  to  the  increase  in 
potential  energy.  As  a  result,  the  total  energy  equation  becomes 


cdT  a  dQ 


(B-7) 


for  an  incompressible  liquid  flowing  vertically  in  a  tube  of  constant 
diameter . 

Assuming  that  no  phase  changes  occur,  an  approximate  energy 
balance  over  the  differential  element  of  depth,  dZ,  yields:  heat  lost 
by  the  liquid  -  heat  transferred  to  the  casing,  or 


dq  -  -WcdT^  -  2irr1U(T1  -  T2)dZ  . 


(B-8) 


The  rate  of  heat  conduction  from  the  casing  to  the  surrounding 
formation  may  be  expressed  as 


dq 


2wk(T2  -  Tc)dZ 


f(t) 


(B-9) 


Eq .  B-9  implies  the  assumption  that  heat  transfers  radially  away  from 
the  wellbore.  The  time  function  f(t)  depends  on  the  conditions 
specified  for  heat  conduction  and  will  be  discussed  later.  Assuming 


the  geothermal  temperature  is  a  linear  function^  of  depth, 


12 


It  is  not  necessary  that  geothermal  temperature  be  linear  with 
depth.  Solutions  may  also  be  obtained  if  geothermal  temperature  is 
represented  graphically  as  a  function  of  depth. 


Te  -  aZ  +  b 


( B - LU ) 


Eqs .  B-9  and  B-10  can  be  substituted  into  Eq .  B-8  to  produce 


+  -  (aZ  +  b)  -  0,  A  i*  0 

&Z  A  A 


(B- 11) 


Wc [k  +  riUf(t) 

A  * 

2jrr^Uk 


(B-12) 


An  integrating  factor  for  Eq.  B-ll  is  eZ/A.  Thus, 


T^eZ/A  -  f(aZ  *  b)eZ/A  dZ  +  C(t)  . (B- 

J  A 


Tiez/A  -  (AZ  -  aA  +  b) eZ//A  +  C(t)  .  (B-14) 


T1(Z,t)  -  aZ  -  aA  +b  +  C(t)e'z/A 


(B- 15) 


The  function  C(t)  can  be  evaluated  from  the  condition  that  T^  -  TQ  for 


C(t)  -  T0( t)  +  aA  -  b 


(B-16) 


And  the  final  expression  for  liquid  temperature  as  a  function  of  depth 
and  time  is 


T^Z.t)  -  aZ  +  b  -  aA  +  [ T0 ( t )  +  aA  -  b]e'z/A 


(B- 17 ) 


To  apply  Eq.  B-17,  it  is  necessary  to  evaluate  the  time  function, 
f(t).  Eq.  B-9  can  be  rearranged  to 


f(t)  -  2*k<T2  *  Tc> 

dq/dZ 


(B- 9a) 


which  is  the  definition  of  this  time  function.  In  this  form,  it  is 
clear  that  the  function  f(t)  has  the  same  relationship  to  transient 
heat  flow  from  a  wellbore  that  the  van  Everdingen-Hurst^  constant 
flux  P(t)  function  has  to  transient  fluid  flow.  In  the  case  of  the 
general  wellbore  heat  transfer  problem,  neither  heat  flux  nor 
temperature  at  the  wellbore  remains  constant  except  in  special  cases. 

A  semi  -  rigorous  treatment  of  transient  heat  conduction  would  involve  a 
complex  superposition  at  each  depth.  Thus,  we  wish  to  find 
approximate  values  of  f(t)  which  will  provide  engineering  accuracy. 
Success  will  be  determined  by  comparison  of  calculated  temperatures 
with  measured  field  temperatures. 


►* 


Fortunately,  many  solutions®  to  transient  heat  and  fluid  flow 
problems  exist  which  may  be  used  to  estimate  f(t).  For  example,  the 
Moss  and  White-*'11  wellbore  heat  transmission  solution  assumes  that 
transient  heat  conduction  to  the  earth  can  be  represented  by  a  line 
source  losing  heat  at  constant  flux.  Carslaw  and  Jaeger-^  present 
graphical  and  analytical  solutions  for  the  cases  of  internal 
cylindrical  sources  losing  heat  at  constant  flux,  constant 
temperature,  and  under  the  radiation  boundary  condition.  Fig.  B-2 
presents  Ramey's  time  function  for  several  different  internal  boundary 
conditions.  As  can  be  seen  from  Fig.  B-2,  the  solutions  presented 
converge  at  long  tines  (about  one  week  or  more) .  This  is  completely 
analogous  to  pressure  build-up  theory  that  at  sufficiently  long  times 
pressure  is  controlled  by  formation  conditions.®  For  times  less 

than  a  dimensionless  time  of  1,000  (i.e.,  at/r2' ^  -  1000),  the 
radiation  boundary  condition  has  been  found  to  produce  logical  values 
for  f(t).  The  radiation  inner  boundary  conditions  is 


6  r 


r-r- 


U2(Ti  -  T2) 


( B - 18) 


where  U2  -  r^U/r2' .  This  boundary  condition  is  comparable  to  the  van 
Everdingen^  skin  effect,  also  well  known  in  pressure  build-up  theory. 
Physically,  Eq .  B-18  states  that  heat  flow  in  the  annular  region 
between  r^  and  r2'  is  controlled  by  steady-state  convection,  rather 


than  conduction. 


conduc 


The  solution  for  this  case  is  presented  by  Carslaw  and  Jaeger^ 
and  is  reproduced  on  Fig.  B-2.  The  time  function  is  seen  to  depend 
upon  (r^U/k) .  However,  the  radiation  boundary  case  does  not  depend 
strongly  upon  (r^U/k)  and  the  solution  to  this  case  approaches  that  of 
the  constant- temperature  cylindrical  source  as  (r^U/k)  approaches 
infinity.  Thus,  the  constant- temperature  cylindrical-source  solution 
is  the  recommended  solution  if  thermal  resistance  in  the  wellbore  is 
negligible.  For  times  greater  than  those  shown  on  Fig.  B-2,  the  line 
source  solution  as  given  by  Eq .  23  is  recommended. 

ASSOCIATED  HEAT  PROBLEMS 

The  solution  presented  by  Eqs .  25  and  26  also  applies  to  wellbore 
heat  problems  other  than  injection  down  tubing.  For  example, 
injection  down  casing  may  be  handled  by  computing  the  over-all 
coefficient  including  only  the  film  coefficient  at  the  casing  wall. 
Wellbore  temperatures  in  a  flowing  well  may  be  computed  if  the  depth 
scale  is  referenced  to  the  producing  interval.  Thus,  T0(t)  becomes 
the  producing  formation  temperature,  and  geothermal  temperature  should 
be  expressed  as  a  function  of  distance  above  the  producing  interval. 

Other  wellbore  heat  problems  can  be  solved  approximately  by 
methods  similar  to  those  used  for  Eq .  25.  That  is,  write  heat 
balances  over  each  flowing  stream  in  the  wellbore  and  assume  that  heat 
loss  from  the  wellbore  may  be  represented  by  Eq .  B-9.  If  two  or  more 
flowing  streams  are  involved,  the  result  will  be  a  higher-order 
differential  equation  than  Eq .  B-ll.  Temperatures  in  each  stream  may 
be  determined,  if  desired.  Note  that  Eqs.  B-8,  B-9,  and  B-10  could 


have  been  solved  for  T2,  the  casing  temperature.  This  problem  may 
have  significance  in  interpretation  of  temperatures  measured  in  the 
annulus  when  fluid  is  flowing  in  the  tubing. 


APPENDIX  C 


RHEOLOGICAL  PRESSURE  LOSS  MODEL  CALCULATIONS 

LAMINAR  FLOW  IN  DRILL  PIPE  AND  ANNULI 

The  drilling  engineer  deals  primarily  with  the  flow  of  drilling 
fluids  down  the  circular  wellbore  of  the  drill  string  and  up  the 
circular  annular  space  between  the  drill  string  and  the  casing  of  the 

Q  C 

open  holeJ  .  If  the  pumping  rate  of  the  drilling  rig's  mud  pumps  is 
low  enough  for  the  flow  to  be  laminar,  the  Newtonian,  Bingham  Plastic, 
or  power  law  model  can  be  employed  to  develop  the  mathematical 
relation  between  flow  rate  and  frictional  pressure  drop.  In  this 
development,  the  following  simplifying  assumptions  are  made^^: 

1.  the  drill  string  is  placed  concentrically  in  the  casing 
or  open  hole , 

2.  the  drill  string  is  not  being  rotated, 

3.  sections  of  open  hole  are  of  circular  shape  and  of  known 

diameter , 

4.  the  drilling  fluid  is  incompressible,  and 

5.  the  flow  is  isothermal. 

In  actuality,  none  of  the  above  simplifying  assumptions  are 
completely  valid^^,  and  the  resulting  system  of  equations  will  not 
perfectly  describe  the  laminar  flow  of  drilling  fluids  in  the  well. 
Also  note  that  the  three  rheological  models  listed  at  the  outset  do 
not  take  into  account  the  thixotropic  nature  of  drilling  mud  and  only 
approximate  the  actual  laminar  flow  fluid  behavior. 


d  «'  v.  •'  -T iTl  s  -*  «T-  -Tj.  -  -  -«  »  *  -  ■*  «  " - 


TURBULENT  FLOW  IN  DRILL  PIPE  AND  ANNULI 


In  many  drilling  operations,  the  drilling  fluid  is  pumped  at  a 
rate  which  is  too  high  for  laminar  flow  to  be  maintained.  The  fluid 
laminae  become  unstable  and  break  into  a  chaotic  diffused  flow 

O  O 

pattern^.  The  transfer  of  momentum  caused  by  this  disorganized  fluid 
movement  causes  the  velocity  distribution  to  become  more  uniform 
across  the  center  portion  of  the  pipe  than  for  laminar  flow.  However, 
a  thin  boundary  layer  of  fluid  near  the  pipe  walls  generally  remains 
in  laminar  flow. 

A  mathematical  development  of  flow  equations  for  turbulent  flow 

O  Q 

has  not  been  possible  to  dateJ  .  However,  a  large  amount  of 
experimental  work  has  been  done  in  straight  sections  of  circular  pipe, 
and  the  factors  influencing  the  onset  of  turbulent  flow  and  the 
frictional  pressure  losses  due  to  turbulent  flow  have  been  identified. 
By  applying  the  method  of  dimensional  analysis,  these  factors  have 
been  grouped  so  that  the  empirical  data  could  be  expressed  in  terms  of 
dimensionless  numbers.  The  onset  of  turbulence  is  associated  with  the 
existence  of  random  fluctuations  in  the  fluid,  and  at  least  on  a  small 
scale,  the  flow  is  inherently  unsteady^.  The  existence  of  turbulent 
flow  can  be  advantageous  in  the  sense  of  providing  increased  heat  and 
mass  transfer  rates.  However,  the  motion  is  extremely  complicated  and 
difficult  to  describe  theoretically. 

For  engineering  purposes,  flow  of  a  Newtonian  fluid  in  pipes 
usually  is  considered  to  be  laminar  if  the  Reynolds  number  is  less 


than  2,100  and  turbulent  if  the  Reynolds  number  is  greater  than  2,100 
However,  for  Reynolds  numbers  of  about  2,000  to  4,000,  the  flow  of 
most  drilling  fluids  is  literally  in  a  transition  region  between 
laminar  flow  and  fully  developed  turbulent  flow. 


SUMMARY  OF  FRICTIONAL  PRESSURE  LOSS  EQUATIONS38 


Newtonian  Fluid  Model 


Mean  Velocity,  v 


v  -  q/2.448d^ 


Annulus 


Flow  Behavior 
Parameters 

Turbulence 

Criteria 


v  -  q/2.448(d2^  -  d^) 


»  “  *300 


NRec  -  2,100 
NRe  -  928pvd/^ 


Annulus 


NRec  -  2,100 

NRe  “  757pv(d2  -  &i)/n 


Laminar  Flow 
Frictional 
Pressure  Loss 


1500d" 


Annulus 

dpf  „  MV 

dL  1000 (d2  -  dx)2 

Pipe 

dp£  ~  p0.75v1.75M0.25 

dL  ISOOd1-25 

Annulus 

dpf  -  p0.75-1.75M0.25 

dL  1 3 96 ( d2  -  d^1-25 

Bingham  Plastic  Model 


Turbulent  Flow 
Frictional 
Pressure  Loss 


Mean  Velocity,  v 


Flow  Behavior 
Parameters 


Turbulence 

Criteria 


Pipe 

v  -  q/2.448d2 
Annulus 


v  -  q/2 . 448 (d2* 


dl2> 


^p  ■  s600  -  hoo 


Pipe 


NHe  -  37,100pryd2/Mp2 


NRec  from  Fig.  4.3338 
NRe  ”  928pvd/Mp 


Annulus 


NHe  -  24,700pry(d2  -  dx)Z/M 


NRe  “  757pv(d2  -  d1)/fip 


Laminar  Flow 
Frictional 
Pressure  Loss 


Turbulent  Flow 
Frictional 
Pressure  Loss 


Power  -  Law  Model 


Mean  Velocity,  v 


Pipe 

dPf  _ 

^PV  + 

ry 

dL. 

1500d2 

225d 

Annulus 

dPf  _ 

4pv 

+ 

ry 

dL 

1000(d2  - 

dx)2 

200(d2 

Pipe 

dPf  -  P0. 75-1. 75^0.25 

dL  1800d125 

Annulus 

dpf  -  p0. 75-1. 75^0. 25 

dL  1396(d2  -  d^1-25 


Pipe 

v  -  q/2 . 448d2 

Annulus 

v  -  q/2 . 448d2 


Flow  Behavior 
Parameters 


n  -  3.32  log(06Oo/03OO) 


Turbulence 

Criteria 


Laminar  Flow 
Frictional 
Pressure  Loss 


Turbulent  Flow 
Frictional 
pressure  Loss 


Nr6C  from  Fig.  4.34^ 


N 


Re 


89,100  p (v) 2 ”n 
K 


0 . 0416d 
.  3+1/n 


Annulus 


N 


Re 


109,000  p(v)2'n 
K 


0.0208(d2  - 
2  +  1/n 


Pipe 


dpf 


Kv11 


'3  +  1/N 
.  0.0416, 


144,000  d1+n 


Annulus 

Kvn 

dpf  -  _ 


3  +  1/N 
.  0.0416, 


144,000  d1+n 


Pipe 

— 9 

dp£  f pwL 

dL~  25. 8d 


Annulus 


dpf 


fpv2 

21 . l(d2  -  dx) 


di) 


dL 


APPENDIX  D 


OVER-ALL  HEAT  TRANSFER  COEFFICIENT 

The  steady- state  rate  of  heat  flow  through  a  wellbore  Q  Btu/hour 
is  proportional  to  the  temperature  difference  between  the  fluid  and 
the  formation,  and  the  cross-sectional  area  perpendicular  to  the 
direction  of  heat  flow^® ,  This  factor  of  proportionality,  which  is 
called  the  over-all  heat  transfer  coefficient,  represents  the  net 
resistance  of  the  flowing  fluid,  tubing,  casing  annulus,  casing  wall, 
and  cement  sheath  to  the  flow  of  heat.  Hence,  we  can  write 

Q  -  UjAjATj  .  (D- 1) 

Eq.  D-l  defines  U,  the  over-all  heat  transfer  coefficient  based  on  the 
characteristic  area  A  and  a  characteristic  temperature  difference  AT. 
The  subscript  j  in  Eq.  D-l  defines  the  surface  area  upon  which  these 
quantities  are  based.  In  theory,  any  radial  surface  could  be  used  to 
determine  the  characteristic  area.  Some  choices  are  more  convenient 
to  work  with  than  others'^.  For  example,  if  a  hot  fluid  were  to  be 
inject  down  tubing,  it  is  preferred  to  let  Aj  be  the  outside  surface 
area  of  an  incremental  length  of  injection  tubing,  2»rrt.0AL,  and  let 
ATj  be  the  difference  between  the  temperature  of  the  flowing  fluid  T^ 
and  the  temperature  at  the  cement- formation  interface  T2 .  The  Uj  - 
UQ ,  referring  to  the  outside  tubing  surface  area,  and  Eq .  D-l  would  be 


Q  -  2ffr  toUCo(Tl  -  T2)AL 


(D-2 


If  the  fluid  is  injected  down  the  casing  or  casing  annulus,  the 
characteristic  area  would  be  the  inside  surface  area  of  the  casing, 
and  Eq .  D-l  would  be  written  as 

Q  -  2,rrciUci(T1  -  T2)AL .  (D-3 

Subscript  ci  refers  to  the  inside  surface  area  of  the  casing. 

An  expression  for  the  over-all  heat  transfer  coefficient  for  any 
well  completion  can  be  found  by  considering  the  heat  transfer 
mechanisms  between  the  flowing  fluid  and  the  cement- formation 
interface.  The  reader  is  referred  to  the  paper  by  Wilhite28  or  any 
good  heat  transfer  textbook  for  a  complete  discussion  and  derivation 
of  the  over-all  heat  transfer  coefficient. 


APPENDIX  E 


PROGRAM  LISTING 


20  COLOR  13,1,5 
30  REM 

********************************************************************* 

40  REM  * 

* 

50  REM  *  HEAT. BAS 

* 

60  REM  *  Heat  Transfer/Temperature  Gradient  Simulation 

* 

70  REM  *  Robert  D.  Pierce 

* 

80  REM  * 

* 

90  REM  *  02  NOVEMBER  1987 

* 

100  REM  * 

* 

110  REM 

********************************************************************* 
120  CLS 

130  KEY  OFF  :  CLEAR 

140  LOCATE  6,12  :  PRINT  CHR$(218) ; 

150  FOR  I  -  1  TO  50  :  PRINT  CHR$(196);  :  NEXT  I 
160  PRINT  CHR$ (191) 

170  FOR  I  -  6  TO  17  :  LOCATE  1+1,12  :  PRINT  CHR$(179)  :  LOCATE  1+1,63 
:  PRINT  CHR$ ( 179 )  :  NEXT  I 

180  LOCATE  8,18  :  PRINT  "A  PC  Simulation  of  Heat  Transfer  and" 

190  LOCATE  10,22  :  PRINT  "Temperature  Distribution  in" 

200  LOCATE  12,25  :  PRINT  "a  Circulating  Wellbore" 

210  LOCATE  14,34  :  PRINT  "by" 

220  LOCATE  16,28  :  PRINT  "Robert  D.  Pierce" 

230  LOCATE  18,12  :  PRINT  CHR$(192); 

240  FOR  I  -  1  TO  50  :  PRINT  CHR$(196);  :  NEXT  I 
250  PRINT  CHR$(217) 

260  PRINT  :  PRINT  :  PRINT 

270  PRINT  TAB(20);"  Press  any  key  to  continue" 

280  IK$-INKEY$  :  IF  IK$-""  THEN  280  ELSE  290 
290  PI  -  3.1415927# 

300  INJ-75  :  REM  ***  PERIOD  OF  INJECTION,  DAYS  *** 

310  KC-25  :  REM  ***  THERMAL  CONDUCTIVITY  OF  CASING,  BTU/HR*FT*DEG  F 

*** 

320  KM— 1 !  :  REM  ***  THERMAL  CONDUCTIVITY  OF  MUD,  BTU/HR*FT*DEG  F  *** 
330  TMI-75!  :  REM  ***  INITIAL  INLET  MUD  TEMP,  DEG  F  *** 

340  NOZ1  -  13  :  REM  ***  NOZZEL  DIAMETER,  32D  OF  INCH  *** 

350  NOZ2  -  13  :  REM  ***  NOZZEL  DIAMETER,  32D  OF  INCH  *** 


360  N0Z3  -  13  :  REM  ***  NOZZEL  DIAMETER,  32D  OF  INCH  *** 

370  D POD-6 . 625  :  REM  ***  DRILL  PIPE  O.D.,  INCHES  *** 

380  DPID-5 . 965  :  REM  ***  DRILL  PIPE  I.D.,  INCHES  *** 

390  DBIT-8 .375  :  REM  ***  DRILL  BIT  O.D.,  INCHES  *** 

400  CSGOD-IO . 7 5  :  REM  ***  CASING  O.D.,  INCHES  *** 

410  CSGWT-51 !  :  REM  ***  CASING  WEIGHT,  LB/FT  *** 

420  CSCID-10!  :  REM  ***  CASING  I.D.,  INCHES  *** 

430  CSGGRD$-"N- 80"  :  REM  ***  CASING  GRADE  *** 

440  Z-15000  :  REM  ***  WELL  DEPTH,  FEET  *** 

450  DELZ-500  :  REM  ***  INCREMENTAL  DEPTH,  FEET  *** 

460  RHOMUD-IO!  :  REM  ***  DENSITY  OF  MUD,  LB/GAL  *** 

470  FR-210  :  REM  ***  FLOW  RATE,  GAL/MIN  *** 

480  MU-110!  :  REM  ***  VISCOSITY,  LB/FT. HR  (CP  -  MU*0.4134)  *** 

490  THETA-70  :  REM  ***  DIAL  READING  ON  ROTATIONAL  VISCOMETER  *** 

500  CPM-.4  :  REM  ***  MUD  SPECIFIC  HEAT,  BTU/LB*DEG  F  *** 

510  CPF- . 2  :  REM  ***  FORMATION  SPECIFIC  HEAT,  BTU/LB . DEG  F  *** 

520  TSA-70  :  REM  ***  SURFACE  AIR  TEMPERATURE,  DEG  F  *** 

530  UP  -  30  :  REM  ***  HEAT  TRANS.  COEFF.  (PIPE),  BTU/HR*SQ . FT*DEG  F 
★  ★★ 

540  UA  -  1!  :  REM  ***  HEAT  TRANS.  COEFF.  (ANNULUS),  BTU/HR*SQ . FT*DEG 
F  *** 

550  GRAD-. 0127  :  REM  ***  ASSUMED  TEMPERATURE  GRADIENT,  DEG  F/FT  *** 

560  TSE  -  59.5  :  REM  ***  SURFACE  EARTH  TEMPERATURE,  DEG  F  *** 

570  SCREEN  0  :  WIDTH  80 
580  COLOR  14,1,1 

590  CLS 

600  GOSUB  920 

610  M$— STRING$ ( 77 , "x" )  :  N$-STRING$ (21 , "x" )  :  R$-STRINC$ ( 26 , "x" ) 

620  PRINT  M$ 

630  PRINT  N$;TAB(57) ;N$ 

640  PRINT  N$ ;TAB( 29) ; "Program  Master  Menu" ;TAB(57) ;N$ 

650  PRINT  N$;TAB(57) ;N$ 

660  PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ;N$ 

670  PRINT  N$;TAB( 57) ;N$ 

680  PRINT  R$ ; "  DATE$  ;  "  xxx  "  ;  TIME$  ";R$ 

690  PRINT 

700  COLOR  4,15  :  PRINT  "(A)  QUIT  :  COLOR  14,1 

710  PRINT  TAB(12);"  (B) .  To  Change/View  Program  Parameters" 

720  PRINT 

730  PRINT  TAB(12);"  (C) . To  Execute  Graphics  Routine" 

740  PRINT 

750  PRINT  TAB(12)  ;  "  (D) . To  Print  Out  Table" 

760  PRINT 

770  PRINT  TAB(  12)  (E) . To  Run  Simulation" 

780  PRINT 

790  PRINT  TAB(12);"  (F) . To  Determine  Pressure  Losses" 

800  PRINT 

810  PRINT  TAB(12);  "To  continue,  press  a  key." 

820  BEEP  :  BEEP 

830  IK$-INKEY$  :  IF  IK$-"”  THEN  GOTO  830  ELSE  GOTO  840 


840 

IF 

IK$-"A" 

OR 

IK$-"a" 

THEN 

END 

850 

IF 

IK$-"B" 

OR 

IK$-"b" 

THEN 

GOTO 

3290 

860 

IF 

IK$-"C" 

OR 

IK$-"c" 

THEN 

GOTO 

6180 

870 

IF 

IK$-"D" 

OR 

IK$-"d" 

THEN 

GOSUB 

.  1780 

880 

IF 

IK$— " E” 

OR 

IK$-"e " 

THEN 

GOTO 

2500 

890 

IF 

IK$-"F" 

OR 

IK$— " f " 

THEN 

GOTO 

4830 

900 

910 

920 

GOTO  570 

REM 

REM  ************  CALCULATIONS  SUBROUTINE 

★  ★★★★★★★★★■*★★★★  •k'k'kic'k'k'k'klr'kic 


930  REM 

990  AREABIT  -  (PI/4096)*( (N0Z1)*2  +  (NOZ2)~2  +  (NOZ3)'2) 

950  PRESBIT  -  ( ( 8 . 311E-05 ) * (RHOMUD)* ( FR) * 2 ) /( ( . 95* . 95 )* ( AREABIT) A 2 ) 
960  IF  Z  <  3000  THEN  KF  -  .6  :  REM  ***  FORMATION  CONDUCTIVITY, 
BTU/HR*FT*DEG  F  *** 

970  IF  Z  >  3000  AND  Z  <-  13000  THEN  KF  -  1.5  +  3 !*(Z/13000) "3 
980  IF  Z  >  13000  THEN  KF  -  4.5 

990  IF  Z  <10000  THEN  RHOFOR  -  122  +  . 0046*Z  :  REM  ***  DENSITY  OF 
FORMATION,  LB/CU.FT  *** 

1000  IF  Z  >  10000  THEN  RHOFOR  -  165 
1010  B  -  ( ( DBIT/2 ) *UA) /( ( DPOD/2 )*UP) 

1020  A  -  ( ( FR*RHOMUD*60) *CPM) /( 2*PI* (DPOD/24 ) *UP) 

1030  Cl  -  (B/(2*A))*( 1+SQR( 1+4/B) ) 

1040  C2  -  (B/(2*A))*(1- SQR( 1+4/B) ) 

1050  C3  -  l  +  ((B/2)*(l+(l+4/B)\5)) 

1060  C4  -  l+( (B/2)*(l- (1+4/B)a . 5)) 

1070  K2  - 


(CRAD*A- (TMI-TSE+-CRAD*A)*EXP(C1*Z)*(1-C3))/(EXP(C2*Z)*(1-C4)-EXP(C1*Z 
)*( 1 -C3) ) 

K1  -  TMI - K2 - TSE+GRAD*A 
OPEN  "0" , 1 , "HEAT1 . DTA" 

FOR  L  -  0  TO  Z  STEP  DELZ 

TP  -  K1*EXP(C1*L)+K2*EXP(C2*L)+GRAD*L+TSE-GRAD*A 
TA  -  K1*C3*EXP(C1*L) +K2*C4*EXP (C2*L) +GRAD*L+TSE 
WRITE  #1,L,TP,TA 
NEXT  L 
CLOSE  1 


1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 

FUNCTION 

1270 


REM 

REM  *** 
REM 

HP-4.36 
OPEN  "0 
FOR  L  - 


BEGIN  RAMEY'S  CALCULATIONS  *** 


, 2 , "RAMEY. DTA" 

0  TO  Z  STEP  DELZ 
UR  -  l/(l/HP+(DPOD-DPID)/( 24*KC ) ) 

FX  -  (  ,04*INJ*24)/((DPOD/2)/12)"2 
FY  -  . 4342945*LOG ( FX ) 

FT  -  ,0018302*FYa3  -  ,045016*FYA2  +  ,49045*FY  -  .44056 
REM  ***  FT  IS  RAMEY'S  TRANSIENT  HEAT  -  CONDUCTION  TIME 
*** 

W  -  1440*FR*RHOMUD 
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1280  R1  -  DPID/2 

1290  REM  ***  KF  -  THERMAL  CONDUCTIVITY  OF  THE  EARTH, 

BTU/DAY*FT*DEG  F 

1300  FIRST  -  W/(2*PI*R1) 

1310  SECOND  -  (CPM*12 ) /( 24*KF*UR) 

1320  THIRD  -  (KF*24)+(RI/2)*UR*FT 

1330  A  -  FIRST*SECOND*THIRD 

1340  T  -  (GRAD*L)+TSA- (GRAD*A)+(TMI+(GRAD*A) -TSA)*EXP( -L/A) 

1350  WRITE  #2 , L, T 

1360  NEXT  L 

1370  CLOSE  2 

1380  RETURN 

1390  REM 

1400  REM 

1410  REM  ***  NEWTONIAN  FLUID  MODEL  *** 

1420  F$  -  "Newtonian" 

1430  VBARDP  -  FR/(2 ,448*(DPID)A2) 

1440  VBARAN  -  FR/ ( 2 . 448*( (CSGIDA2 ) - ( DPID' 2 ) ) ) 

1450  NREDP  -  928*RHOMUD*VBARDP*(DPID) /(MU* .4134) 

1460  NREAN  -  757*RHOMUD*VBARAN*(CSGOD- DPID) /(MU* . 4134) 

1470  PR  -  CPM*( ,4134*MU)/KM 

1480  IF  NREDP  <  2100  THEN  FLOCON$  -  "LAMINAR" 

1490  IF  NREAN  <  2100  THEN  CONFLO$  -  "LAMINAR" 

1500  IF  NREDP  >-  2100  THEN  FLOCON$  -  "TURBULENT" 

1510  IF  NREAN  >-  2100  THEN  CONFLO$  -  "TURBULENT" 

1520  IF  NREDP  <  2100  THEN  1550  ELSE  1640 
1530  IF  NREAN  <  2100  THEN  1550  ELSE  1640 
1540  REM 

1550  REM .  LAMINAR  FLOW  CALCULATIONS 


1560  FRIC  -  64/NREDP 

1570  NUBAR  -  4.36 

1580  HP  -  (NUBAR*KM) /DPID 

1590  DPRESP  -  ((MU*.4134)*VBARDP*Z)/(1500*(DPID'2)) 

1600  DPRESA  -  ( (MU* ,4134)*VBARAN*Z)/(1000*( (CSGID-DPOD) *2) ) 
1610  DPTOT  -  DPRESP  +  DPRESA  +  PRESBIT 
1620  GOSUB  920 
1630  GOTO  5120 


1640  REM . TURBULENT  FLOW  CALCULATIONS 


1650  REM 

1660  IF  NREDP  <  20000  THEN  FRIC  -  . 316*NREDP' ( - . 25)  ELSE  GOTO  1670 
1670  FRIC  -  ,184*NREDP'(- .2) 

1680  NUBAR-  . 023*(NREDP' . 8 ) *PRa ( 1/3) 


1690  HP  -  (NUBAR*KM)/DPID 
1700  DPRESP  - 

((RHOMUDA . 75 ) * (VBARDP A 1 ,75)*((MU*.4134)A.25) *Z) /( 1800*(DPIDA 1.25)) 
1710  DPRESA  - 

( (RHOMUD' . 75) * (VBARAN A 1 . 75 ) * ( (MU* . 41 34) A . 25 )*Z)/( 1396* (DPID' 1 . 25) ) 
1720  DPTOT  -  DPRESP  +  DPRESA  +  PRESBIT 


1730  GOSUB  920 
1740  GOTO  5120 
1750  REM 

1760  REM  ***********  END  CALCULATIONS 
**************************************** 

1770  REM 

1780  REM  ***********  ROUTINE  TO  PRINT  FILE  CONTENTS 
************************* 

1790  REM 
1800  CLS 

1810  INPUT  "Do  you  want  the  Holmes  &  Swift  model  (Y/N)" ;Y$ 

1820  IF  LEFT$(Y$, 1)-"Y"  OR  LEFT$ (Y$ , l)-"y"  THEN  1830  ELSE  2180 
1830  PRINT  :  PRINT 

1840  PRINT  :  INPUT  "DO  YOU  WANT  A  HARD  COPY  (Y/N)”;Y$ 

1850  IF  LEFT$(Y$,1)-"Y"  OR  LEFT$(Y$ , l)-"y"  THEN  2030 
1860  PRINT  :  PRINT 

18  70  A$  -  "###.##'•  ;  B$  -  "######.##" 

1880  PRINT  TAB(3) ; "DEPTH" ;TAB(23) ; "DRILL  PIPE  TEMPERATURE";  TAB(50); 
"ANNULUS  TEMPERATURE" 

1890  PRINT  TAB(3)  ;  " . ";TAB(23);" . "  ;  TAB(50); 


1900  OPEN  " I " , 1 , "HEATl . DTA" 

1910  IF  EOF ( 1 )  THEN  CLOSE  1  :  GOTO  1990 

1920  INPUT  #1 , L , TP , TA 

1930  PRINT  USING  "#####  FEET" ; L; 

1940  PRINT  TAB (26) ; "  "; 

1950  PRINT  USING  "###.##  DEG  F" ; TP ; 

1960  PRINT  TAB(51) ; "  "; 

1970  PRINT  USING  "###.##  DEG  F" ; TA 
1980  GOTO  1910 

1990  PRINT  :  PRINT  "HIT  ANY  KEY  TO  RETURN  TO  MAIN  MENU" 

2000  IK$  -  INKEY$  :  IF  IK$  -  ""  THEN  GOTO  2000  ELSE  GOTO  2010 
2010  IK$  -  INKEY$  :  IF  IK$  -  "INKEY$"  THEN  GOTO  570 
2020  RETURN 
2030  REM 

2040  REM  ***********  ROUTINE  TO  PRINT  A  HARD  COPY  OF  TABLE 
★★★★★★★★★★★★★★★★★ 

2050  REM 

2060  LPRINT  :  LPRINT 

2070  LPRINT  TAB(3) ; "DEPTH" ;TAB(23) ; "DRILL  PIPE  TEMPERATURE";  TAB(50) 
"ANNULUS  TEMPERATURE" 


2080 

M 

LPRINT  TAB( 3 ) ; " . 

tl 

" ; TAB( 2  3 ) ; " 

2090 

OPEN  "I”,l, "HEATl. DTA" 

2100 

2110 

IF  EOF(l)  THEN  CLOSE 
INPUT  #1 , L, TP , TA 

1  : 

RETURN 

2120 

2130 

LPRINT  USING  "###*# 
LPRINT  TAB (26) ; "  " ; 

FEET" 

;  L; 

2140 

2150 

LPRINT  USING  "###.## 
LPRINT  TAB (51) ; "  "  ; 

DEG 

F" ; TP ; 
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2160  IF  LEFT$ (Y$ , 1)-"Y"  OR  LEFT? (Y? , l)-ny"  THEN  2540  ELSE  570 
2170  GOTO  2100 

2180  REM  ***  ROUTINE  TO  PRINT  RAMEY'S  DATA  *** 

2190  CLS 

2200  PRINT  :  PRINT 

2210  PRINT  :  INPUT  "DO  YOU  WANT  A  HARD  COPY  (Y/N)";Y$ 

2220  IF  LEFT? (Y?,1)-"Y"  OR  LEFT? (Y? , l)-"y"  THEN  2370 
2230  PRINT  :  PRINT 

2240  PRINT  TAB(3) ; "DEPTH" ;TAB(23) ; "WELLBORE  FLUID  TEMPERATURE" 

2250  PRINT  TAB  ( 3 )  ;  " . ";TAB(23);" . " 

2260  OPEN  "I" , 1, "RAMEY. DTA" 

2270  IF  EOF ( 1)  THEN  CLOSE  1  :  GOTO  2330 

2280  INPUT  #1 ,L,T 

2290  PRINT  USING  "#####  FEET" ; L; 

2300  PRINT  TAB (26) ; "  " ; 

2310  PRINT  USING  "###.##  DEG  F";T 
2320  GOTO  2270 

2330  PRINT  :  PRINT  "HIT  ANY  KEY  TO  RETURN  TO  MAIN  MENU" 

2340  IK?  -  INKEY?  :  IF  IK?  -  ""  THEN  GOTO  2340  ELSE  GOTO  2350 
2350  IK?  -  INKEY?  :  IF  IK?  -  "INKEY?"  THEN  GOTO  570 
2360  RETURN 
2370  REM 

2380  REM  ***********  ROUTINE  TO  PRINT  A  HARD  COPY  OF  TABLE 
***************** 

2390  REM 

2400  LPRINT  :  LPRINT 

2410  LPRINT  TAB(3) ; "DEPTH" ;TAB(23) ; "WELLBORE  FLUID  TEMPERATURE" 

2420  LPRINT  TAB ( 3 )  ;  " . ";TAB(23);" . " 

2430  OPEN  "I" ,1, "RAMEY. DTA" 

2440  IF  EOF ( 1 )  THEN  CLOSE  1  :  RETURN 

2450  INPUT  #1 , L, T 

2460  LPRINT  USING  "#####  FEET" ; L; 

2470  LPRINT  TAB (26) ; 

2480  LPRINT  USING  "###.##  DEG  F";T 
2490  GOTO  2440 

2500  REM  **************  SIMULATION  ******************************** 
2510  CLS 

2520  INPUT  "Have  you  run  the  fluid  model  yet  (Y/N)";Y? 

2530  IF  LEFT?(Y? , 1)-"Y"  OR  LEFT?(Y? , l)-"y"  THEN  2540  ELSE  4830 


2540  SCREEN  1,0  :  COLOR  14,0 
2550  GOSUB  2620  :  REM  ---  DRAW 
2560  LOCATE  4,2  :  PRINT  "Press 
2570  GOSUB  3150  :  REM  ---  DRAW  ’ 
2580  IK?  -  INKEY?  :  IF  IK?  -  "" 
2590  GOSUB  2650  :  REM  ---  UPDAT 
2600  GOSUB  2870  :  REM  ---  SIMUL 
2610  GOTO  570 

2620  REM  ---  DRAW  BORDER  --- 
2630  LINE  (1,1)  -  (318,198)  ,2, B 
2640  RETURN 


COLOR  14,0 

REM  ---  DRAW  BORDER  --- 
PRINT  "Press  any  key  to  continue." 

REM  ---  DRAW  WELLBORE  --- 

:  IF  IK?  -  ""  THEN  GOTO  2580  ELSE  2590 

REM  ---  UPDATE  SCREEN  --- 

REM  ---  SIMULATION  EXECUTION  --- 
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2650  REM  ---  PRINT  DATA  ON  SCREEN  --- 

2660  LOCATE  7,2  :  PRINT  "FLOW  RATE  ; FR ; "  GPM" 

2670  LOCATE  8,2  :  PRINT  "VBARDP  :  PRINT  USING  "####.»#  " ; VBARDP ; : 
PRINT  "  FPS" 

2680  IF  FLOCON$  -  "TURBULENT"  THEN  R$  =  "TURB"  ELSE  R$-"LAM" 

2690  LOCATE  10,2  :  PRINT  "VBARAN  =" ; :  PRINT  USING  "###«.*#  " ; VBARAN ; : 
PRINT  "  FPS" 

2700  IF  CONFLO$  -  "TURBULENT"  THEN  Q$  -  "TURB"  ELSE  Q$-"LAM” 

2710  LOCATE  9,2  :  PRINT  "NREDP  :  PRINT  USING 
" ; NREDP; : PRINT  R$ 

2720  LOCATE  11,2  ;  PRINT  "NREAN  — " ; :  PRINT  USING  "###,###,### 

" ; NREAN; : PRINT  Q$ 

2730  LOCATE  12,2  :  PRINT  "DRILL  PIPE  TEMP  -  ";;PRINT  USING  "###.w 
"; TP;; PRINT  "DEG  F" 

2740  LOCATE  13,2  :  PRINT  "ANNULUS  TEMP  -  ";: PRINT  USING  "*##.# 

";TA;; PRINT  "DEG  F” 

2750  LOCATE  17,2  :  PRINT  "PRESS;" 

2760  LOCATE  19,3  :  PRINT  " (S)  CONTINUE  SIMULATION" 

2770  LOCATE  20,3  :  PRINT  "(A)" 

2780  LOCATE  21,3  :  PRINT  "(B)" 

2790  LOCATE  22,3  ;  PRINT  " (M)  RETURN  TO  MAIN  MENU" 

2800  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  2800  ELSE  GOTO  2810 

2810  IF  IK$-"S"  OR  IK$-"s"  THEN  GOTO  2860 

2820  IF  IK$-"A"  OR  IK$-"a"  THEN  GOTO  2860 

2830  IF  IK$-"B"  OR  IK$-"b"  THEN  GOTO  2860 

2840  IF  IK$-"C"  OR  IK$-"c"  THEN  GOTO  2860 

2850  IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  570 

2860  RETURN 

2870  REM  ---  SIMULATION  EXECUTION  --- 

2880  X  -  232 

2890  FOR  Y  -  20  TO  180 

2900  LINE  (X- 7 , Y)  -  (X+7 , Y+10) , 2 , BF 

2910  LINE  (X-7.Y-10)  -  (X+7,Y),0,BF 

2920  FOR  N  -  1  TO  100 

2930  NEXT  N 

2940  NEXT  Y 

2950  XI  -  220  :  X2  -  244 

2960  FOR  Y  -  189  TO  100  STEP  -1 

2970  LINE  (X1-3.Y)  -  (Xl+3 , Y-10) , 2 , BF 

2980  LINE  (XI - 3 , Y+10)  -  (Xl  +  3 , Y) , 0 , BF 

2990  LINE  (X2 - 3 , Y)  -  (X2+3 , Y- 10) , 2 , BF 

3000  LINE  (X2 -  3 , Y+10)  -  (X2  +  3 , Y) , 0 , BF 

3010  FOR  N  -  1  TO  100 

3020  NEXT  N 

3030  NEXT  Y 

3040  Xl-215  :  X2-249 

3050  FOR  Y  -  100  TO  19  STEP  -1 

3060  LINE  (X1-6.Y)  -  fXl+6 , Y - 10) , 2 , BF 

3070  LINE  (XI  -  6 , Y+10)  -  (Xl+6 , Y) , 0 , BF 

3080  LINE  (X2-6.Y)  -  (X2+6 , Y- 10) , 2 , BF 


3090  LINE  (X2-6.Y+10)  -  (X2+6 , Y) ,  0 ,  BF 
3100  FOR  N  -  I  TO  100 
3110  NEXT  N 

3120  NEXT  Y 

3130  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  3130  ELSE  GOTO  3140 
3140  RETURN 

3150  REM  ---  DRAW  WELLBORE  --- 
3160  LINE  (208,10)  -  (208, 100), 2 
3170  LINE  -  (216,100) ,2 
3180  LINE  -  (216, 190), 2 
3190  LINE  -  (248, 190), 2 
3200  LINE  -  (248,100) ,2 
3210  LINE  -  (256, 100), 2 
3220  LINE  -  (256,10) ,2 
3230  LINE  (224,10)  -  (224, 180), 1 
3240  LINE  (240,10)  -  (240, 180), 1 
3250  LINE  (208,10)  -  (224, 10), 2 
3260  LINE  (240,10)  -  (256, 10), 2 
3270  LINE  (200,0)  -  (200,199) 

3280  RETURN 

3290  REM  ***  PARAMETER  CHANGE  MENU  *** 

3300  CLS 
3310  PRINT  M$ 

3320  PRINT  N$ ; TAB (57) ; N$ 

3330  PRINT  N$ ;TAB(26) ; "Parameter  Change/View  Menus” ; TAB(57) ;N$ 

3340  PRINT  N$; TAB (57) ;N$ 

3350  PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ,N$ 

3360  PRINT  N$; TAB (57) ;N$ 

3370  PRINT  R$ ; "  DATE$  ;  "  xxx  "  ;  TIME$  ; "  ";R$ 

3380  PRINT 

3390  COLOR  4,15  :  PRINT  "(A)  QUIT  :  COLOR  14,1 

3400  PRINT  TAB (12) ; "  (B) . To  Change  Drilling  Fluid 

Parameters" 

3410  PRINT 

3420  PRINT  TAB (12) ; "  (G) . To  Change  Drill  Pipe  Parameters 

3430  PRINT 

3440  PRINT  TAB(12);"  (D) . To  Change  Wellbore/Formation 

Parameters" 

3450  PRINT 

3460  PRINT  TAB(12);"  (E) . To  Change  Casing  Parameters" 

3470  PRINT 

3480  PRINT  TAB (12) ; "  (M) . To  Return  To  Main  Menu" 

3490  PRINT  :  PRINT  TAB(12);  "To  continue,  press  a  key."  :  BEEP 

3500  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  3500  ELSE  GOTO  3510 

3510  IF  IK$-"A"  OR  IK$-"a"  THEN  END 

3520  IF  IK$-"B"  OR  IK$-"b"  THEN  GOTO  3570 

3530  IF  IK$-"C"  OR  IK$-"c"  THEN  GOTO  3910 

3540  IF  IK$-"D"  OR  IK$-"d"  THEN  GOTO  4200 

3550  IF  IK$-"E"  OR  IK$-"e"  THEN  GOTO  4510 

3560  IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  570 


3570  REM  ***  DRILLING  FLUID  PARAMETERS  *** 

3580  CLS 
3590  PRINT  M$ 

3600  PRINT  N$ ;TAB(57) ;N$ 

3610  PRINT  N$;TAB(27) ; "Drilling  Fluid  Parameters" ; TAB (57) ;N$ 

3620  PRINT  N$;TAB(57) ;N$ 

3630  PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ;N$ 

3640  PRINT  N$; TAB (57) ;N$ 

3650  PRINT  R$ ; "  "  ;  DATE$  ;  "  xxx  "  ;  TIME$  ;  "  " ;R$ 

3660  PRINT 

3670  COLOR  4,15  :  PRINT  "(A)  QUIT  :  COLOR  14,1 
3680  PRINT  TAB (12) ; "  (B)  Thermal  Conductivity  of  Mud  - 

" ;KM; "Btu/(hr*f t*deg  F) " 

3690  PRINT 

3700  PRINT  TAB (12) ; "  (C)  Initial  Inlet  Mud  Temperature  -  ";TMI;"deg 

F" 

3710  PRINT 

3720  PRINT  TAB (12) ; "  (D)  Density  of  Mud  -  " ; RHOMUD ; "  Lb/gal" 

3730  PRINT 

3740  PRINT  TAB (12) ; "  (E)  Flow  Rate  -  "  ;  FR ; "  Gal/min" 

3750  PRINT 

3760  PRINT  TAB (12) ; "  (F)  Absolute  Viscosity  of  Mud  -  " ;MU;" 

Lb/(ft*hr) " ; "  ("MU*. 4134"  CP  )" 

3770  PRINT 

3780  PRINT  TAB (12) ; "  (G)  Specific  heat  of  Mud  -  ”;CPM;"  Btu/(lb*deg 

F) " 

3790  PRINT 

3800  PRINT  TAB(12);"  (M)  Return  to  Parameters  Menu" 

3810  PRINT  :  PRINT  TAB(12);  "To  continue,  press  a  key."  :  BEEP 
3820  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  3820  ELSE  GOTO  3830 
3830  IF  IK$— "A"  OR  IK$-"a"  THEN  END 

3840  IF  IK$-"B"  OR  IK$-"b"  THEN  INPUT  "New  value  for  thermal 
conductivity  of  mud  -";KM  :  GOTO  3570 

3850  IF  IK$-"C"  OR  IK$-"c"  THEN  INPUT  "New  value  for  initial  inlet 
mud  temperature  -";TMI  :  GOTO  3570 

3860  IF  IK$-"D"  OR  IK$-"d"  THEN  INPUT  "New  value  for  mud  density  = 
RHOMUD  :  GOTO  3570 

3870  IF  IK$-"E"  OR  IK$-"e"  THEN  INPUT  "New  value  for  flow  rate  -  ” ; FR 
:  GOTO  3570 

3880  IF  IK$-"F"  OR  IK$-"f"  THEN  INPUT  "New  value  for  viscosity  -  ";MU 
:  GOTO  3570 

3890  IF  IK$— "G"  OR  IK$-"g"  THEN  INPUT  "New  value  for  specific  heat  of 
mud  -  " ; CPM  :  GOTO  3570 

3900  IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  3290 
3910  REM  ***  DRILL  PIPE  PARAMETERS  *** 

3920  CLS 
3930  PRINT  M$ 

3940  PRINT  N$;TAB(57) ;N$ 

3950  PRINT  N$;TAB(29) ; "Drill  Pipe  Parameters" ;TAB(57) ;N$ 

3960  PRINT  N$;TAB( 57) ;N$ 
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3970  PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ;N$ 

3980  PRINT  N$;TAB(57) ;N$ 

3990  PRINT  R$ ; "  "  ;  DATE$  ;  "  xxx  "  ;  TIME$  ;  "  ";R$ 

4000  PRINT 

4010  COLOR  4,15  :  PRINT  "(A)  QUIT  ” ;  :  COLOR  14,1 

4020  PRINT  TAB(12) ; ”  (B)  Heat  Trans.  Coeff.  at  Pipe  Wall  - 

” ;UP; "Btu/(hr*ftA2*deg  F)" 

4030  PRINT 

4040  PRINT  TAB(12);"  (C)  Drill  Pipe  Inner  Diameter  -  DPID inches" 

4050  PRINT 

4060  PRINT  TAB(12);"  (D)  Drill  Pipe  Outer  Diameter  -  "; DPOD; "inches" 


(E)  Incremental  Length  -  ”;DELZ;"  ft" 


(M)  Return  to  Parameter  Menu" 

"To  continue,  press  a  key. 


BEEP 


coefficient 
4160  IF  IK$ 
diameter  ; DPID 
4170  IF  IK$-"D" 
diameter  =-";DPOD 


4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 


New  value  for  heat  transfer 
New  value  for  drill  pipe  inner 
New  value  for  drill  pipe  outer 
New  value  for  delz  -  ” ; DELZ  : 


4070  PRINT 

4080  PRINT  TAB (12) ; ” 

4090  PRINT 

4100  PRINT  TAB (12) ; " 

4110  PRINT  :  PRINT  TAB(12); 

4120  PRINT  :  PRINT 

4130  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  4130  ELSE  GOTO  4140 
4140  IF  IK$— "A"  OR  IK$-"a"  THEN  END 

4150  IF  IK$-"B"  OR  IK$-"b"  THEN  INPUT 

" ;UP  :  GOTO  3910 
"C"  OR  IK$-"c"  THEN  INPUT 
:  GOTO  3910 
OR  IK$-"d"  THEN  INPUT 
:  GOTO  3910 

4180  IF  IK$-"E"  OR  IK$-"e"  THEN  INPUT 
GOTO  3910 

IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  3290 
REM  ***  WELLBORE/ FORMAT ION  PARAMETERS  *** 

CLS 

PRINT  M$ 

PRINT  N$;TAB(57) ;N$ 

PRINT  N$;TAB(29) ; "Wellbore/Formation" ; TAB( 57 ) ;N$ 
PRINT  N$;TAB(33) ; "Parameters" ;TAB(57) ;N$ 

PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ;N$ 

PRINT  N$;TAB(57) ;N$ 

PRINT  R$ ; "  "  ;  DATE$  ;  "  xxx  "  ;  TIME$  ;  "  ";R$ 
PRINT 

COLOR  4,15  :  PRINT  "(A)  QUIT  ” ;  :  COLOR  14,1 


iriod  of  Injection  -  " ; INJ ; "  days" 
tmperature  Gradient  -  ”;GRAD;"  deg  F/ft" 
irface  Air  Temperature  -  ";TSA;"  deg  F" 
illbore  Diameter  -  ";DBIT;"  inches" 
iat  transfer  coefficient  (annulus)  -";UA;" 

(H)  Formation  Density  -  ";RHOFOR;"  Ib/cu  ft" 


4310 

PRINT  TAB (12)  ;  " 

(B) 

4320 

PRINT  TAB (12) ; " 

(C) 

4330 

PRINT  TAB (12) ; " 

(D) 

4340 

PRINT  TAB (12) ; " 

(E) 

4350 

PRINT  TAB (12)  ;  " 

(F) 

4360 

PRINT  TAB (12) ; " 

(G) 

Btu/(hr*sq  ft*deg  F)" 

4370 

PRINT  TAB (12) ; " 

(H) 
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4380  PRINT  TAB (12) ; "  (I)  Formation  Conductivity  -  " ; KF ; " 

Btu/(hr*f t*deg  F) " 

4390  PRINT  TAB(12);"  (M)  Return  to  Parameter  Menu" 

4400  PRINT  :  PRINT  TAB(12);  "To  continue,  press  a  key.”  :  BEEP 
4410  PRINT  :  PRINT 

4420  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  4420  ELSE  GOTO  4430 

4430  IF  IK$-"A"  OR  IK$-"a"  THEN  END 

4440  IF  IK$-"B"  OR  IK$-"b"  THEN  INPUT  "New  value  for  well  depth  -";Z 

:  GOTO  4200 

4450  IF  IK$-"C"  OR  IK$-"c"  THEN  INPUT  "New  value  for  period  of 
injection  -  ” ; INJ  :  GOTO  4200 

4460  IF  IK$-"D"  OR  IK$-”d"  THEN  INPUT  "New  value  for  temp  grad  - 
" ; GRAD  :  GOTO  4200 

4470  IF  IK$-"E"  OR  IK$-"e"  THEN  INPUT  "New  value  for  surface  air  temp 
-  " ; TSA  :  GOTO  4200 

4480  IF  IK$-"F"  OR  IK$-"f"  THEN  INPUT  "New  value  for  wellbore 
diameter  —  ";DBIT  :  GOTO  4200 

4490  IF  IK$-"G"  OR  IK$-"g"  THEN  INPUT  "New  value  for  heat  transfer 

coeff.  (annulus  -  ";UA  :  GOTO  4200 

4500  IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  3290 

4510  REM  ***  CASING  PROGRAM  PARAMETERS  *** 

4520  CLS 
4530  PRINT  M$ 

4540  PRINT  N$; TAB (57) ;N$ 

4550  PRINT  N$;TAB(27) ; "Casing  Program  Parameters" ;TAB(57) ;N$ 

4560  PRINT  N$; TAB (57) ;N$ 

4570  PRINT  N$;TAB(33) ; "Select  One" ;TAB(57) ;N$ 

4580  PRINT  N$;TAB(57) ;N$ 

4590  PRINT  R$ ; "  "  ;  DATE$  ;  "  xxx  "  ;  TIME$  ;  "  ";R$ 

4600  PRINT 

4610  COLOR  4,15  :  PRINT  "(A)  QUIT  :  COLOR  14,1 

4620  PRINT  TAB(12);"  (B)  Thermal  Conductivity  of  Casing  - 

" ;KC; "Btu/(hr*ft*deg  F)" 

4630  PRINT 

4640  PRINT  TAB(12) ; "  (C)  Inside  Diameter  of  Casing  -  ";CSGID;" 

Inches" 

4650  PRINT 

4660  PRINT  TAB (12) ; "  (D)  Outside  Diameter  of  Casing  -  ";CSGOD;" 

Inches" 

4670  PRINT 

4680  PRINT  TAB (12) ; ”  (E)  Weight  of  Casing  -  ";CSGWT;”  Lb/ft" 

4690  PRINT 

4700  PRINT  TAB (12) ; "  (F)  Grade  of  Casing  -  ";CSGGRD$ 

4710  PRINT 

4720  PRINT  TAB(12);"  (M)  Return  to  Parameter  Menu" 

4730  PRINT  :  PRINT  TAB(12);  "To  continue,  press  a  key."  :  BEEP 
4740  PRINT  :  PRINT 

4750  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  4750  ELSE  GOTO  4760 
4760  IF  IK$-"A”  OR  IK$-"a”  THEN  END 
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4770  IF  IK$-"B"  OR  IK$-"b"  THEN  INPUT  "New  value  for  thermal 
conductivity  of  casing  -";KC  :  GOTO  4510 

4780  IF  IK$— "C"  OR  IK$-"c"  THEN  INPUT  "New  value  for  casing  I.D.  - 
” ;CSGID  :  GOTO  4510 

4790  IF  IK$-"D"  OR  IK$-"d"  THEN  INPUT  "New  value  for  casing  O.D.  - 
" ; CSGOD  :  GOTO  4510 

4800  IF  IK$-"E"  OR  IK$-"e"  THEN  INPUT  "New  value  for  casing  weight  - 
" ; CSGWT  :  GOTO  4510 

4810  IF  IK$-"F"  OR  IK$-"f"  THEN  INPUT  "New  value  for  casing  grade  - 
";CSGGRD$  :  GOTO  4510 

4820  IF  IK$— "M"  OR  IK$-"m"  THEN  GOTO  3290 
4830  REM  ***  MUD  RHEOLOGY  MODEL  *** 

4840  CLS 
4850  PRINT  M$ 

4860  PRINT  N$;TAB(57) ;N$ 

4870  PRINT  N$;TAB(27) ; "Fluid  Rheology  Model" ;TAB(57) ;N$ 

4880  PRINT  N$;TAB(57) ;N$ 

4890  PRINT  N$ ; TAB (33) ; "Select  One" ;TAB(57) ;N$ 

4900  PRINT  N$;TAB(57) ;N$ 

4910  PRINT  R$ ; "  "  ;  DATE$  ;  "  xxx  "  ;  TIME$  ;  "  ";R$ 

4920  PRINT 

4930  COLOR  4,15  :  PRINT  "(A)  QUIT  ";  :  COLOR  14,1 
4940  PRINT  TAB (20) ; "  (B)  Newtonian" 

4950  PRINT 

4960  PRINT  TAB(20) ; "  (C)  Bingham  Plastic" 

4970  PRINT 

4980  PRINT  TAB (20) ; "  (D)  Power  -  Law" 

4990  PRINT 

5000  PRINT  TAB (20) ; "  (E)  Return  to  simulation" 

5010  PRINT 

5020  PRINT  TAB (20) ; "  (M)  Return  to  Main  Menu" 

5030  PRINT  :  PRINT  TAB(12);  "To  continue,  press  any  key."  :  BEEP 
5040  PRINT  :  PRINT 

5050  IK$-INKEY$  :  IF  IK$-""  THEN  GOTO  5050  ELSE  GOTO  5060 

5060  IF  IK$-"A"  OR  IK$-"a"  THEN  END 

5070  IF  IK$-"B"  OR  IK$-"b"  THEN  GOTO  1410 

5080  IF  IK$-"C"  OR  IK$-"c"  THEN  GOTO  5560 

5090  IF  IK$— "D"  OR  IK$-"d"  THEN  GOTO  5790 

5100  IF  IK$-"E"  OR  IK$— " e "  THEN  GOTO  2500 

5110  IF  IK$-"M"  OR  IK$-"m"  THEN  GOTO  570 

5120  CLS 

5130  PRINT  M$ 

5140  PRINT  N$;TAB(57) ;N$ 

5150  PRINT  N$;TAB(28) ; "Fluid  Flow  Conditions" ; TAB( 57 ); N$ 

5160  PRINT  N$;TAB( 57) ;N$ 

5170  PRINT  R$ ; "  ";  DATE$  ;  "  xxx  "  ;  TIME$  ;"  ";R$ 

5180  PRINT 

5190  COLOR  4,15  :  PRINT  "Fluid  model  chosen  is  " ; F$  :  COLOR  14,1 
5200  PRINT 
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5210  PRINT  TAB( 10) COLOR  14 , 0 : PRINT"Average  velocity  in  drill  pipe  - 
PRINT  USING  "#«#».##  " ; VBARDP ; :  PRINT  "  ft/sec" 

5220  COLOR  14,1,1 

5230  PRINT  TAB ( 10) COLOR  14 , 0 : PRINT"Average  velocity  in  the  annulus 
-  PRINT  USING  ”#**#.#*  " ; VBARAN ; :  PRINT  "  ft/sec” 

5240  COLOR  14,1,1 

5250  PRINT  TAB( 10) ;: COLOR  14 , 0 : PRINT"Reynold' s  Number  in  drill  pipe  - 
PRINT  USING  "###,##*,##*  " ; NREDP ; :  PRINT  FLOCON$ 

5260  COLOR  14,1,1 

5270  PRINT  TAB( 10) COLOR  14 , 0 : PRINT"Reynold' s  Number  in  annulus  - 
PRINT  USING  " ### , ### , " ; NREAN ; :  PRINT  CONFLO$ 

5280  COLOR  14,1,1 

5290  PRINT  TAB( 10) COLOR  14,0: PRINT"Darcy  friction  factor  - 
PRINT  USING  "#.#####  " ; FRIC 
5300  COLOR  14,1,1 

5310  PRINT  TAB( 10) COLOR  14 , 0 : PRINT"Prandtl  Number  -  PRINT  USING 

"##.###  ";PR 

5320  COLOR  14,1,1 

5330  PRINT  TAB( 10) COLOR  14 , 0 : PRINT"Nusselt  Number  -  PRINT  USING 
"##,###.##  " ; NUBAR 
5340  COLOR  14,1,1 

5350  PRINT  TAB ( 10) COLOR  14,0:PRINT  "Pressure  loss  in  drill  pipe  - 
PRINT  USING  "##,###.#  DPRESP PRINT  "  psi" 

5360  COLOR  14,1,1 

5370  PRINT  TAB( 10) COLOR  14,0:PRINT  "Pressure  loss  in  annulus  - 
■' ;  :  PRINT  USING  "#,###.#  DPRESA PRINT  ”  psi" 

5380  COLOR  14,1,1 

5390  PRINT  TAB( 10) COLOR  14,0:PRINT  "Pressure  loss  in  bit  -  PRINT 
USING  "#,###.#  " ; PRESBIT ; : PRINT  "  psi" 

5400  COLOR  14,1,1 

5410  PRINT  TAB( 10) COLOR  14,0:PRINT  "Total  Pressure  loss  in  system  — 
" ; : PRINT  USING  "##,###.#"; DPTOT PRINT  "  psi" 

5420  COLOR  14,1,1 

5430  IF  IK$-"C"  OR  IK$-"c"  THEN  PRINT  TAB ( 10) COLOR  14,0: PRINT 
"Critical  Reynold's  number  (drill  pipe)  -  " ;NRECDP  ELSE  5460 
5440  COLOR  14,1,1 

5450  PRINT  TAB( 10) ;: COLOR  14,0: PRINT  "Critical  Reynold's  number 

(annulus)  -  " ; NRECAN  :  COLOR  14,1,1 

5460  IF  IK$— "D"  OR  IK$-”d"  THEN  GOTO  5470  ELSE  5510 

5470  PRINT  TAB( 10) ;: COLOR  14,0: PRINT  "Flow  behavior  index  (n)  - 

"  ;  NLAW 

5480  COLOR  14,1,1 

5490  PRINT  TAB( 10) ;: COLOR  14,0:PRINT  "Consistency  index  (K)  -  " ;KLAW 
5500  COLOR  14,1,1 

5510  PRINT  :  PRINT  "Hit  any  key  to  return  to  Fluid  Rheology  menu." 
5520  BEEP 

5530  IK$-INKEY$  :  IF  IK$  -  ""  THEN  GOTO  5530  ELSE  GOTO  5540 
5540  GOTO  4830 
5550  REM 

5560  REM  ***  BINGHAM  PLASTIC  MODEL  *** 
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5570  REM 

5580  F$  -  "Bingham  Plastic" 

5590  VBARDP  -  FR/(2 ,448*(DPID) A2) 

5600  VBARAN  -  FR/<2 ,448*( (CSGIDA2) - (DPIDA2) ) ) 

5610  MUP  -  THETA  -  (MU*. 4134) 

5620  YP  -  (MU*. 4134)  -  MUP 

5630  NHEDP  -  ( 37100 ! *RHOMUD*YP*(DPIDA2) )/(MUPA2) 

5640  NHEAN  -  (24700*RHOMUD*YP*( (CSGID-DPOD) A2) )/(MUPA2) 

5650  NREDP  -  (928*RHOMUD*VBARDP*DPID)/MUP 

5660  NREAN  -  ( 757*RH0MUD*VBARAN*(CSGID-DP0D) )/MUP 

5670  HED1  -  LOG (NHEDP) 

5680  HED2  -  LOG (NHEAN) 

5690  NRECDP  -  EXP( - . 0019841*(HED1A3)  +  . 080586*(HED1A2)  - 
(,71107*HED1)  +  9.4408) 

5700  NRECAN  -  EXP( - . 0019841*(HED2 A3)  +  . 080586*(HED2A2)  - 
( ,71107*HED2)  +  9.4408) 

5710  PR  -  CPM*( . 4134*MU) /KM 

5720  IF  NREDP  <  NRECDP  THEN  FLOCON$  -  "LAMINAR" 

5730  IF  NREAN  <  NRECAN  THEN  CONFLO$  -  "LAMINAR" 

5740  IF  NREDP  >-  NRECDP  THEN  FLOCON$  -  "TURBULENT" 

5750  IF  NREAN  >-  NRECDP  THEN  CONFLO$  -  "TURBULENT" 

5760  IF  NREDP  <  NRECDP  THEN  1550  ELSE  1640 
5770  IF  NREAN  <  NRECAN  THEN  1550  ELSE  1640 
5780  REM 

5790  REM  ***  POWER  -  LAW  MODEL  *** 

5800  REM 

5810  F$  -  "Power  Law" 

5820  VBARDP  -  FR/( 2 . 448*(DPID) A2 ) 

5830  VBARAN  -  FR/(2 . 448*( (CSGIDA2 ) - (DPIDA2) ) ) 

5840  NLAW  -  3 . 32* . 4342945*LOG(THETA/(MU* . 4134) ) 

5850  KLAW  -  (510*(MU*,4134) )/(511aNLAW) 

5860  NREDP  - 

( (89100 !*RH0MUD*( VBARDP A( 2 -NLAW) ) /KLAW) )*( ( ( . 0416*DPID)/(3+(l/NLAW) ) ) - 
aNLAW) 

5870  NREAN 

-(( 109000 !*RHOMUD*VBARANa(2-NLAW)/KLAW)*( (( . 0208*(CSGID-DPOD) )/(2+(l/- 


NLAW)) 

)aNLAW)) 

5880 

PR  -  CPM*( ,4134*MU)/KM 

5890 

IF  NREDP 

<  2100  THEN  FL0C0N$  - 

"LAMINAR" 

5900 

IF  NREAN 

<  2100  THEN  CONFLO$  - 

"LAMINAR" 

5910 

IF  NREDP 

>-  2100  THEN  FL0C0N$  -  "TURBULENT 

5920 

IF  NREAN 

>-  2100  THEN  CONFLO$  -  "TURBULENT 

5930 

IF  NREDP 

<  2100  THEN  5950  ELSE 

6050 

5940 

IF  NREAN 

<  2100  THEN  5950  ELSE 

6050 

5950 

REM 

5960 

REM - 

----  LAMINAR  FLOW  CALCULATIONS 

5970 

FRIC  -  64/NREDP 

5980 

NUBAR  -  4.36 

5990 

HP  -  (NUBAR*KM) /DPID 

6000  DPRESP  - 

((KLAW*(VBARDPaNLAW)*Z*((3+1/NLAW)/. 0208 )aNLAW)/( 144000 ! *DPIDA ( l+NLAW 

))) 

6010  DPRESA  - 

( (KLAW* (VBARANaNLAW) *Z*( ( 2+1/NLAW) / . 0208 ) aNLAW) /( 144000 ! * (CSGID-DPOD) 
A( l+NLAW))) 

6020  DPTOT  -  DPRESP  +  DPRESA  +  PRESBIT 
6030  GOSUB  920 
6040  GOTO  5120 

6050  REM . TURBULENT  FLOW  CALCULATIONS 


6060  REM 

6070  IF  NREDP  <  20000  THEN  FRIC  -  .  316*NREDPA ( - . 25 )  ELSE  GOTO  6100 
6080  IF  NREAN  <  20000  THEN  FRIC  -  . 316*NREANA ( - . 25)  ELSE  GOTO  6100 
6090  FRIC  -  . 184*NREDPA(- .2) 

6100  FRIC  -  . 184*NREANa(- .2) 

6110  NUBAR  -  ,023*(NREDPa ,8)*PRa(1/3) 

6120  HP  -  (NUBAR+KM) /DPID 

6130  DPRESP  -  (FRIC*RHOMUD*Z*VBARDPA2)/(25.8*DPID) 

6140  DPRESA  -  (FRIC*RHOMUD*Z*VBARANA2)/(21 . l*(CSGID-DPOD) ) 

6150  DPTOT  -  DPRESP  +  DPRESA  +  PRESBIT 

6160  GOSUB  920 

6170  GOTO  5120 

6180  RUN  "QKPLT.EXE" 
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